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SUMMARY 

In this report we consider a liquid in a cylindrical container having 
either circular or rectangular cross-section. The bottom of the container 
is a rigid, perfect heat conductor while the upper surface is a free interface 
between the liquid and the bounding gas. The bottom of the container is 
uniformly heated and heat is lost to the gas at the top. The basic state 
consists of zero velocity (motionless liquid) having a linear, purely con- 
ductive temperature profile. If the vertical temperature gradient is large 
enough, this pure conduction state becomes unstable due to the joint effects 
of two mechanisms. The surface tension on the interface depends on tempera- 
ture and Marangoni convection sets in. The non-dimensional measure of the 
suriace-tension gradient is the Marangoni number M, The vertical gravity 
field sets up an adverse density gradient so that buoyancy effects lead to 
convection. The non-dimensional measure of the buoyancy is the Rayleigh 
number R. 

We formulate the general problem of nonlinear convective instability driven 
by the joint effects of thermocapillarity and buoyancy. The upper free sur- 
face has a general heat transfer condition applied and the interface is 
allowed to deform. Sidewalls confine the layer. The problem solved 
involves a special case of the general problem. The surface tension on the 
interface is so large that surface deflections are neglected. The sidewalls 
are adiabatic and impenetrable but for mathematical simplicity are allowed 
to be "slippery’*. 

In order to determine the effect of sidewalls on the critical 
Marangoni number in finite containers, it is first necessary to develop 


tile basic equations and boundary conditions. The basic equations are the 
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Navier* Stokes , continuity and energy equations^ and the boundarv conditions 
express the conditions of no slip and cot^stancy of temperature at the bottom 
surface of the layer, conservation of momentum and energy at the free inter- 
face, and the sidewall conditions alluded to above. 

If we denote by a the ratio of the mean fluid depth to the radius of 
the cylinder, then we find the critical condition's H > 
instability of tiie pure conduction state. At most vitlues of the aspect 
ratio a, linear stability theory selects a single cell shape of steady 
Marangoni convection that replaces the pure condition when the conduction 
is unstable. The cell shape is characterized by an integer m since all 
solutions are proportional to exp(ini0) where 0 is the cylindrical azimuthal 
angle. Thus m ■« 0 corresponds to axisymmetric ring cells. Alternatively, 
m » 1 corresponds to convection patterns in which there is up flow in one 
half the container and downflow in the other half. Clearly, the modes for 
larger m correspond to more complicated cellular patterns. The linear 
theory shows that as the aspect ratio a ii. increased from zero that the 
preferred modes follow the sequence m* 1, m* 2, m *0, m* 3, etc. We 
study convection near intersections between modes and away from such 
intersections. 

The linear instability theory for the case of the rectangular container 
is even more complicated, as two aspect rotios are necessary to describe 
the container. We denote the dimensionless Length and width of the rectangl 
as (a^,a,,). Again, modes of convection are denoted by integers (m^,m^) 
corresponding to the number of cells in the horizontal directions (x,y). 

Thus (1,0) is a single roll cell with axis in the y direction. Parameter 
Studies result in a map in the (a^,a^) plane which gives the modes which 


-Ic- 


are preferred according to linear theory. We shall be especially concerned 
with boxes which are of shallow depth relative to their length a^; this 
means we shall focus for fixed a^ on the progression of roll cells (m,0) as 
the length a^ of the box increases. Again, we study convection near inter- 
sections between modes and away from these. 

In all the above calculations, both for the circular and the rectangular 
geometries, results are obtained for various Rayleigh numbers (measuring 
buoyancy) and various Biot numbers (measuring heat transfer to the gas). 
However, for all the nonlinear results obtained both the Rayleigh and Biot 
numbers are set to zero. Subsequent work will cover cases involving non- 
zero values of these numbers. 

In order to detexnnine the mode of convection that is observable in 
an experiment, it is necessary to develop a nonlinear theory. Such a theory 
predicts the amplitude and direction of fluid motion, and for container 
sizes for where two different modes of convection compete, is capable of 
making predictions of the convection pattern as a result of such a nonlinear 
competition. 

Consider a single mode, whose amplitude is denoted by A(t). Then 
according to linear theory, A satisfies an equation of the form, 

(3.1) 

where v is a constant which depends upon the mode in question, the aspect 

ratio, and the Prandtl number. According to (3.1), A will grow in time 

v;hen M > M , and decay when M < M : this is the result of linear theory, 
c c 

UTien M (3.1) is no longer valid tor all times, and a nonlinear theory 
is necessary to describe the convection. The nonlinear analog to eq. (3.1) 
may be derived using eigenfunction expansions. In the s:implest case, 


-Id- 


such a procedure results in the antpllcude equation. 


^ (M-M )A - kA^ 


(3.2) 


where k is a computable constant which depends upon the mode, the aspect 
ratio, and the Prandtl number* Conclusions regarding the form and ampli- 
tude of the convection may be made by examining the solutions to eq. (3.2) 
and their stability. For example, there are three solutions to eq. (3.2), 

A ^ 0, A » ± [ k"^(M - M^)] It is easy to show that for k > 0, M^, the 

null solution is unstable and the other solutions are stable. Conversely, 

if k < 0, real solutions exist for M < M only, and these are all unstable. 

c 

Thus for k > 0, the system evolves to a stare of steady convection whose 
magnitude is given by [aI. The purpose of our nonlinear theory is to derive 
the amplitude equations analogous to eq. (3.2) for Marangoni convecti n in 
finite containers and to determine the steady solutions and their stability. 
For example, (subcritical) convection can in some cases exist for M< M , 
even though a linear theory would predict no convection. 

When the container has aspect ratio close to those values tor which 
two modes become unstable simultaneously, one must write a pair of equa- 
tions for the amplitudes, A^(t) and A.,(t), say. Tliese will be of the form 


dA 




(3.3a) 


dA, 
^2 dt 




(3.3b) 


where f^, f, are nonlinear functions of their arguments, end f^(0,0) = 
f^(0,0) ^ 0. Clciv’ly, near these aspect ratios, more steady solutions 
exist to eq. (3.3^ then in the case of a simple mod', and the ^'‘ehavior is 
more complex. 


k 


le 


We have derived the forms of eqs, (3.3) for finite amplitude Marangoni 
convection in cylinders and boxes, and have determined completely the steady 
solutions and their stability. We discuss in detail here the results for 
the cylindrical containers. Results for the rectangular box are similar 
in many ways. 

A nonlinear theory of convection is developed so that at each aspect 
ratio a of the cylinder, we can compute the whole fluid and temperature 
field. As long as m 0 so the convection is non-axisyrametric , we find 
that the instability is supercritical i.e. even when disturbances of con- 
duction become larger, there is no sustained convection for M < M^(a). 

However, when m » 0 and the convection is axisymmetr ic , we find a range of 
subcritical convection at M < long as the disturbance level is 

large enough. The theory for m = 0 also distinguishes the flow direction at the 
cylinder center. It rises in the center when the Prandtl number Pr of 
the fluid satisfies Pr ^ L and descends in the center of Pr is small. 

There is a special value of the aspect ratio a where two distinct 

instability modes, m = 1 and m * 2, of linear theory are equally likely. 

In this case our nonlinear theory can be applied for, say, a near and 
predictions can be made for fixed a and for M increasing above M . We 

find very interesting phenomena* For a >a , the linearized prediction (at 

A 

M “ M ) of convection is mode m ■ 2. Then, as M is raised, there can be 
c 

a sudden transition to m » 1 convection and possibly time-periodic motion. 
Alternatively, the mode tn * 2 can persist. On the other hand for a < a^, 
the linearized prediction (at M “ M^) of convection is mode m * 1 and the 
nonlinear theory shows chat mode m = 1 must be replaced by mod;^ m 2. We 
thus see a lack of synmotry in the behavior on the two sides oi a “ a^. 
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There is a special value of a^ of the aspect ratio a where two distinct 

instability modes, m ■ 2 and m ■ 0, of linear theory are equally likely. 

For a > a , linearized theory gives axisv-mraetric convection as preferred. 

Our nonlinear theory shows that this convection is subcritical (i,e, occurs 

at H < As M is increased, there is a tendency to remain on the 

axisymnetric mode with no further transition though m * 2 convection can 

exist under certain conditions* On the other hand, when a < a^, there is 

a very complicated sequence of transitions predicted* At the neutral curve 

(M * M^) linearized theory states that mode m * 2 appears* However, nonlinear 

theory gives the following result. As M is Increased above 

mode m = 2 convection begins. As M is increased further, this mode m « 2 

becomes unstable and there is a sudden transition to the axisywmetric mode 

m = 0 which then persists for increasing M* However, if M now 

decreases below M there would be a sudden transition not to the inode 
c 

m » 2 but to pure conduction* Hence, there can be dynamic hysteresis 
loops in this case. Notice again the lack of syrmetry for a >a^ and a < a . 

The behaviors are quite different* 

These behaviors are predicted to occur for Marangoni convection, but 
not for buoyancy-driven convection* In all our nonlinear work we have neglected 
buoyancy by setting the Rayleigh number R » 0 and have examined pure 
Marangoni convection M 0* Rosenblat, in an independent study, has 
examined the complementary problem of pure buoyancy driven convection, 

R f 0, M * 0, in a cylinder. He finds linear stability curves similar to those 
obtained here. He analyzes the nonlinear theory and finds the following 
behaviors. All modes of convection, including the mode m » 0 , bifurcate 
supercrit ically . Near intersections, the transitions for R increasing 
are always symmetric ir. the sense that on either side, the mode predicted 
by linear theory (at R - R^/) suddenly become unstabu (at R " R to mixed 
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aodes composed of the modes at the Intersection point. Hence, there Is no 
subcrltical instability, no lack of synnetry near intersections and no 
dynamical hysteresis. 

Comparison of the above two sets of predictions is crucial for 
evaluation of possible space experiments since under micro-gravity condi- 
tions the first behaviors should occur, while for thick layers on Earth, 
buoyancy dominates and the second behaviors should occur. 

The structure of this report is as follows. Chapter I gives the full 
analysis for convection in circular cylinders* This includes discussion of 
previous work, formulation of the governing equations and boundary conditions, 
linear and nonlinear stability theory, and finally discussion of the physical 
implications of the theory. Chapter II gives the parallel development for 
convection in rectangular cylinders. Certain mathematical details are 
postponed to Appendices. Thus, all Tasks required by the contract are 
ftilfilled though in a different order than specified. The structure was 
chosen in order to emphasize the physical understanding of the phenomena. 

All aspects of the circular geometry are complated first. Then all aspects 
of the rectangular geometry are examined. 


CHAPTER I . CONVECTION W CIRCULAR CYLINDERS 

I, INTROIXJCTION 


Consider a uniform layer of liquid having infinite horizontal extent, 
bounded on the bottom by a solid plate and having a free surface on the top. 
When the plate is heated with respect to the gas at the free surface, a 
purely conductive static state may exist, in which P is the (constant) magni- 
tude of the temperature gradient. If the free surface possesses surface 
tension a, the variations a(T) of surface tension with temperature T can induce 
Marangoni instability . This thermocapillary instability was identified and 
explained by Pearson (1958) who showed, using a linear stability theory, .that 
a critical value of Marangoni number M must be exceeded before the conductive 
state becomes unstable. Here, 


M 


do 

dT 


Pd^ 


iC 


cf 0 


( 1 . 1 ) 


where — is the (negative) rate of change of surface tension with temperature, 
the subscript zero denoting a constant value at a given reference temperature 
Tq, d is the thickness of the undisturbed layer; and are the thermal 
diffusivity and dynamic viscosity of the liquid, respectively. 

The critjxal value of M depends on other parameters: the surface Biot 

number h, which is the non-dimensional version of the heat transfer coefficient 
at the free surface, and a capillary number C, which is a non-dimensional 
version of the mean surface tension ^ aCF^). Here 




( 1 . 2 ) 


Pearson (1958) limited his analysis to C ^ 0, which means that the free surface 
does not deform as a result of disturbances. He finds, for h = 0 and a perfectly 
conducting lower boundary, ^ 79.6 in the absence of gravity. 
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When a vertical gravity field is present, so taat the layer is heated 
from below, Nield (1964) shows, using Pearson's model, that buoyancy-induced 
instabilities and thermocapillary instabilities reinforce on another. 

Since the work of Pearson and Nield, linear stability theory on the 
Marangoni convection problem has been extended in several directions to include 
two dynaraical phases, C 0, and further effects of an imposed vertical gravity 
field; Scriven and Sternling (1964), Smith (1966), Zeren and Reynolds (1972). 

A recent survey, Sorensen (1978), discusses these as well as many other 
extensions. Palmer and Berg (1971) find that the theory of Nield (1964) pre- 
dicts well the experimental conditions for the onset of convection in shallow 
layers. 

As in tiie case of Ray leigh-Benard convection due to buoyancy effects, the 

horizontal planform of the convective state and the amplitude of the motion 

(and heat transfer) are undetermined by linear stability theory. In addition 

as in Rayleigh-Benard convection, the critical point M , corresponding to the 

c 

critical horizontal wave number, is infinitely degenerate; there are an 
infinite number of planforms allowable by linear theory. Nonlinear effects 
presumably select from this set those that appear in experiments. The first 
nonlinear analysis of Marangoni Instability is due to Scanlon and Segei (1967), 

They consider an inf inite-Prandti-number fluid, an infinity-deep layer and 
only twc ^ planform functions from the infinite set. Their prediction is that 
hexagonal-cells is the only planform that exists aiid is stable when the 
conductive state becomes unstable. Hexagons exist and are stable for an interval 
of M M so that subcriticai convection is oredicted. They do not attemnc to 
enlarge the set of planform functions beyond the chosen two. However, their 
prediction is in qualitative accord with experimental observations, Koschmieder (1967) 
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always finds very regular hexagons in shallow layers » the regularity stenming 
from the precise thermal controls used* These pictures should be compared 
with those of Benard (1900) who finds irregular, polygonal cells in convection 
dominated by thermocapillarity* 

The only other nonlinear analysis of Marangoni instability is due to 
Kraska and Sani (1979)* They consider six planform functions (Including those 
of Scanlon and Segel) and also complete a nonlinear analysis. However, they 
encounter difficulties in analyzing the stability of their nonlinear states 
and find a lack of closure in adding any seventh state to the original six* 

These difficulties may be related Co the adjoint operator they define, which 
seems not to be appropriate to the problem. 

All of the above work is stimulated by the recognition that Marangoni 
instability is the prototype instability associated with the transport of 
heat and/or mass across interfaces,. The presence of such instabilities can 
augment the transport rates by orders of magnitude compared to that of pure 
conduction. Particular ajpects of these trar. sport processes have motivated 
many of the generalizations mentioned above. 

In recent years there has been a new interest in Hows driven by thermo- 
capillary forces* With the advent of Spacelab research, there is the possibility 
of performing fluid mechanics experiments in space. The micro-gravity environ- 
ment allows interfacial forces to dominate gravity so chat one may be able to 
design experiments that focus on effects submerged on Carth by gravitational 
effects. The present work is motivated by such corisiderations. 

There are special difficulties in inicrogravi ty environments. Those 
are illustrated in Figure 1. 


Figure 1 


Here, a mass of liquid is placed in a wide dish. On Earth, as shown in 
Fig. la, the liquid fills the dish nearly uniformly except for small menisci 
near the edges. A meniscus has a width W that scales as B ^d where B is 
the Bond number, 

B = pQgd^/aQ , (1.3) 

Pq is the liquid density and g is the gravitational potential. Thus on Earth 
W is usually a fraction of d, whereas for gravity at 10 ^g, the meniscus ’’fills 
the whole dish”, as shown in Fig. lb. Alternatively, depending upon the 
liquid volume, the contact angle and the wetting properties of the solid, the 
liquid mass might form in to a sessile drop or sphere or be detached entirely 
from the container as shown in Figures Ic and Id. Thus, any Marangoni insta- 
bility study to be applied to microgravity conditions must allow for the 
experimental necessity of near-by, confining sidewalls. Figure 2a shows a 
possible configuration of the free surface wbon the aspect ratio is near 
unity. The sidewalls help to both confiiu^ the liquid and limit the degree of 
curvature of the fre* surface. 


Figure 2 


- 6 - 







In the present studies we address the problem of Marangoni instabilities 
in a cylinder. Our particular interest is to explore the nonlinear i *ter- 
actions near M = and in particular transitions from one convecti.\/e state 
to another that may occur at supercritical conditions. Given the difficulty 
of such a study, we make several simplifying assumptions, (i) We let the 
capillary number C 0. H^nce, the top free surface is non-deformable. In 
addition, we take the contact angle at the sidewalls to be compatible with 
a flat free surface. Thus, in the basic, conductive state there are no 
menisci and in the convective state the free surface remains flat (Figure 2b). 

( 11 ) In order to allow the linear stability theory to be solved using normal 
moaes (separation of variables), we idealize the sidewall boundary conditions 
in the following way. The sidewalls consist of a circular cylinder through 
which there is zero heat flux and zero mass flow. However, we allow the 
walls to be ‘'slippery** so that rather than the no~slip condition, we apply 
the conditions of zero tangential vorticity. Clearly, such an idealization 
modifies the predictions of a theory. We shall discuss in detail some 
implications of this idealization and suggest how results of such a theory 
should be applied. 

In the language of bifurcation theory, we study the perfect problem. 
Imperfections due to horizontal temperature gradients or free-surface deflections 
will not be treated here. Note that the change from **slippery" to no-slip 
sidewalls is not an imperfection . 

Given the model described, we shall find the nonlinear Marangoni 
convective-states and determine their stability. Clearly, the presence of 
the sidewall makes the sp ctruin at M = M discrete so that the complete behavior 
can be examiaeu. Of course, hexagonal cells for small containers will never 
apj^ear, since the allowal le cell shapes are luT.naLGd by the sidewall constraints. 

' L PAGE IS 

quality 
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We obrain its behavior* its amplitude and hence can find all transport quantities 
of the convection. At certain aspect ratios two linearized modes are equally 
likely at critical conditions. We analyze such double eigenvalues and find 
certain strong behaviors. The transitions as M is increased for aspect ratios 
on one side of the double eigenvalues, differ substantially from those on 
the other side of the double eigenvalue. Such a demarcation of behaviors 
should be characteristic of the nonlinear ities and hence should be observable 
In an experiment* 

The technique of nonlinear stability theory we use is due to Rosenblat 
(1979), who makes an eigenfunction expansion of the nonlinear problem. This 
’*infinite-matrix*' form is systematically simplified by defining a new sniall 
parameter related to the separation of the eigenvalues of the matrix. The 
results coincide with the usual waekly nonlinear bifurcation theory very near 
M , but gives a wider range of validity. It is the wider range that allows 
us to predict the successive transitions. 
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2, FORMULATION 


Consider a viscous liquid , which partially fills a cylindrical 
container of circular cross-section, llxe mean depth of the liquid 1 
d» and the radius of the cylinder is taken to be ad, so that the mea:; aspect 
ratio (ratio of radius to mean depth) is a. The axis of the cylinder is anti- 
parallel to the direction of gravity, and the upper surface of the liquln is 
open to an ambient gas* 

The liquid is assumed to be Newtonian, to have constant 
viscosity Mq and to be heat conducting with constant thermal 
dlffusivlty The density p* of the liquid is taken to vary with the 

temperature according to the equation of state 

- P^[l - a(T* - T^)] (2.1) 

where p^, Tq are constant reference density and temperature respectively, and 
where a Is the coefficient of volumetric expansion. The liquid-gas interface 
has a surface tension a* which varies with temperature according to the 
formula 


0^ - o, (T* - 
0 I s 


^ 0 ^ 


( 2 . 2 ) 


where a , a are constants and T* Is the temperature at the Interface. 

0 1 ^ 

The tJoussinesq approximation is assumed; the governing equations in 
L ' bulk of the liquid are the Navler-S tokes , conrlnulty and energy-balance 


equations : 




+ (y*‘V*)v*} - 7 **t* - o*gx 


(2.3) 


7**y* m Q 


(2 4) 


4^— + (v*«V*)T* - K- 
3t* - 0 


(2.5) 


where v* Is the velocity vector, g is the acceleration dae to gravity, is 
the stress tensor, z is a unit vector in the upward vertical direction, and 
t*, refer to dlraensici;al time and space derivatives respectively. For a 
Newtonian liquid the stress tensor has the form 


* -p*I ^ + (V*v^) 1 


( 2 . 6 ) 


where p* Is the pressure and I Is the Identity tensor. 


We shall work in a cylindrical polar i:oordinate system, with the 
origin at the center of the lower circular boimdviry. 'Hie mean height of the 
liquid tii located at 2* * d, and the lateral boundary at ^ aa. 

The lower boundary is at z* - 0 and is assumed to be a rigid, 
perfectly conducting plane. Fhus, we have 


r* *» Tq + 3d on “ 0 


10 


and 


V* ■ 0 on z* ■ 0 


(2.7b) 


where we shall take 3 > 0. 

The upper boundary is a free surface, denoted S, which is cooled hy 
heat transfer to the gas. We take the heat transfer relationship to have the 
form 


-n*V*T* - h*(T* - Tq) + 3 on S (2.8a) 

where n is unit normal to S directed into the gas, and h* is the heat 
transfer coefficient. We write 


S: z* » d* + l*(l*» 


(2.9) 


where r^ denotes position vector In the horizontal plane, so that the 
kinematic boundary condition takes the form 


- + v*‘7*(z* - d* - n*) 

o t * 


0 on S. 


(2.8b) 


Finally, we have the dynamic (stress) boundary conditions, which are 


T**n “ 2H*o*a - n * (n * 7*o*) on S 


(2.8c) 


where H* is the mean curvature of the surface. Equation (2.8c) expresses the 


fact that the stress tensor experiences a jump In the normal direction due to 
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surface Cension and a Jump In the tangential direction due to (thermally- 
induced) variations in surface tension. 

The lateral boundary la assumed to be adiabatic. This give the 
condition 


V*T**n - 0 on S 


(2.10a) 


where Is unit nornal vector in the outward direction. If the boundary were 
rigid we would have v* * ^ there- However, as explained in the Introduction, 
we take instead the idealized, mathematically simpler condition that the 
sidewall is a nondef ormable surface on which the tangential vorticity is 
zero. Nondeformability implies that coincides with r* ■ ad, 

0 v' 2* < d + n*, and we have the conditions 


v*-^ » 0 on S 


(2.10b) 


and 


To close the system we require a condition at the contact line where 
the free surface meets the sidewalk- In this paper we shall assume the 
contact angle to be fixed at the value tj/2, so that we have the condition 




X V* X V*) = 0 on S 


(2.10c) 


7*n**n 


0 . 


(2.U) 


kn iqulllbrl'.:m solution of the system (2.1)-(2.1l) Is the follow’ng; 
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V* - 0, T* - T - Tq + 0(d - z*) , 

p* - p - Pq fl - a 0(d - z*)l , 0 * - Oq, .1* - 0, (2.12) 

T* “ T - pggd i(z*/d - 1) + j a^d(z*/d - 

This is the conduction solution, whose stablity we propose to exaiaine. The 
gas pressure has here been set equal to tero* Note that the llquid**ga8 
Interface la flat, and that the temperature at the upper boundary is Tq, which 
implies that 0 is the vertical temperature gradient. 

We now perturb the basic state (2«12), and at the same time introduce 
approp ’Imenslonless forms of the equations and boundary conditions. We 

scale lengths on the depth d, and write 

r = , 2 = z*/d , = *^*/d (2.13) 

with r - (r, z) in cylindrical polar coordinates. Unit vectors in the 
corresponding directions will be denoted (r, z). The liquid occupies the 

region 0<r<a, 0<z<14-nin this dimensionless coordinate system; the 
equation of the free surface S is 

2 - I + n(r^. c) (2.14) 

where r^^ is dimensionless position vector in the horizontal plane and t is 
dimensionless time, defined below; the lateral boundary Is r - a, 


0 < 2 < 1 n. 
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Since our Interest Is focussed on motions driven by surface-tension 
gradients, it is appropriate to base velocities on the Marangoni velocity 
scale V^, defined by 


V - 0 3d/u . (2.15) 

The time-scale will then be d/Vj^^, so that we write 

t - t*V /d. (2.16) 

M 

We now introduce the following representations: 

V* - V V - V (u, V, w) 

T* - T + 3d8 “ Tq + 3d(l - z + 0) 

p* - p - ap|^d6 - p^ - ap^3d(l - z + 9) (2. 1 7) 

0* - + Oj^3d(n - 9) 

T* - T + (u^V^/d)! • P^gd 1(Z - 1) + J a3d(z - l)^}l + (UQV^/d)T, 

where the perturbation field quantititea 9, v and t are functions of r and t, 
and where the forms of p* and o* derive directly from (2.1) and (2.2) 
respectively. Tlie stress tensor t can be written in the form 

i " -pL + (Vv + C7v)^l 


( 2 . 18 ) 
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where p Is the fressure. 

We substitute (2.13)-(2.18) into the governing equations and boundary 
conditions (2.I)-(2.IL). T> e equations In the bulk of the liquid for the 
perturbation field quantities are found to be 

1 2 —1 
Pr M lgY+ (v*7)y| - -Vp + V y + M R0z (2.19) 


?• V ■ 0 


( 2 . 20 ) 


M - w + (yV)0[ - 7 0 


where the parameters appearing here are 


Marangonl number 


Rayleigh number 


Prandtl number 


M 




Oj^3d 


O0gd‘^p 


P - 


"o^o 


Pr 


0 0 


( 2 . 21 ) 


(2.22a) 


(2.22b) 


(2.22c) 


The boundary conditions are as follows. On the lower horizontal 
boundary, equations (2.7) reduce to 


0*0 on z*0, 0<.riia, 


(2.23a) 


V =• 0 on z = 0, 


0 \ 'S 3 . 


(2.23b) 
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On the upper surface S, given now by (2.14), the heat loss condition (2.8a) 
becomes 


where 


n*V0 • n »2 


1 - h(0 - n) on S 


(2.2<*a) 


h =* h”d (2.22d) 

is the surface Biot number. Tlie kinematic condition (2.8b) becomes 


iil + Jh V 

3t 3r r 3(j> 


w 


0 on S . 


(2.2'4b) 


It Is convenient to decompose the stress condition (2.8c) into an equaticn 
normal to the surface and an equation tangential to the surface. If we write 
normal and tangential components as 


(T-n) H (T*n)-n, (T*n) = (T*n) - ((T*n)*n]n, 

as-Q - "-ta " - 

we obtain the dynamic surface conditions 

G(n + I u6Tn^) + MC(T*n)^ - 2H 1 1 + MC( n - 8) } - 0 on S (2.24c) 


and 


(T*n) + n * (n V(n - a)| » 0 on S (2.24d) 

a - ta ■* ” 


where H ■ H^d ts dimensionless mean curvature, and where the two new 

parameters appearing here are ^ 

Bond number C ■ — — (2,22w) 

— 0 ^ 
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caplllary number 


V‘ 


(2.22f) 


On the lateral boundary we have that ■ r. Hence the boundary 
conditions (2.10) become 


ii 

3r 


0 on 


(2.25a) 


u ■ 0 on Sj^ 


(2.25b) 



3w 

3r 


0 on S . 

La 


C.25c) 


The contact-angle condition (2.11) becomes 

■ 0 on r = a. (2.26) 

Dr 

Finally, we recall Che stipulation that the mean depth of the liquid 
is d. This is effectively a condition of volume conservation and can be 
expressed as 


a 2x 

/ J rn(r, 9 , t)d;dr - 0. (2.27) 

0 0 

The problem to be studied comprises the system of equations (2.1*))- 
(2.21), together with the conditions (2.23)-(2.27). There are six parameters 

M, R, Pr, G, h and C, defined by equations (2.2?.); the aspect ratio a is 
a paramcCer of the problv\n,» 


For our purposes the principal par.i.ueter Is the Marangonl number M, 




-17- 




whlch aeasures the effect of surface-tension gradients. We shall determine 
Mg, the critical value of M at which the conduction solution becomes unstable, 
as a function of the other parameters, that is. 




M - M (R, Pr, G, C, h, a). (2.28) 

c c 

We shall then Investigate the properties of the convection as M increases 
above for various values of the other parameters. 
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3. ZERO CAPILLARY NllMBl-R LIMIT 

In this paper wc shall confine our analysis co the case where the 
capillary number C is zero. llic capillary number is associated with 
deflection of the free surface; the limir C ♦ 0 corresponds to a flat 
surface, Davis and Homsy (1980), a restriction which enables onset of con- 
vection to be analyzed relatively simply. In this limit equation (2.24c) 
reduces to 


G(n + I o^dri^) - 211 


(3.1) 


where the standard definition of the mean curvature H is given by 


2H 




(3.2) 


being the gradient In the horizontal plane). Equation (3*1) Is thus a 
differential equation for the surface deflection n - n(r, t), which is 
required to satisfy the conditions (2.26) and (2.27). A solution Is clearly 


n 


H 0. 


(3.1) 


representing an undeformed free surface. The fact that there can be r.o otl.er 
solution In the class of functions n with |nl sufficiently small is easily 
established with the aid of the iipplicit lunction theoren. We confine our 
atlentiv^n to these weakly-nonlinear interactions that ^pply when .M is e Icse 
to and when convective "i.^tions are of small ainplitude. 
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On the basis of this reasoning we infer that the upper free surface 
reaaina flat and undeformed. The equation of S is now z»l,0<r<a, and 
the outward unit normal to It Is n • z. The lateral boundary is r * a, 

0 < z < 1. The problem to be solved in the limit C 0 therefore comprises 
the equations (2. 19)-'(2.2l) In the bulk of the liquid, together with boundary 
conditions as follows. On the lower boundary the conditions ire (2.23) or, 
equivalently. 


3«^«v»w«0 on 2 * 0 , 0<r<a. (3.4) 


On the upper free surface (2.24a) reduces to 


~ hO * 0 on z > 1, 0 < r < a (3.5a) 

3z 


while (2.24b) becomes simply 


w»0 on z«l, 0sr<a. O. >b) 


In view of (3.3) the condition (2.24c) is redundant. We simplify (2.24d) by 


( TMl) 


ta 




^ (rw)]' 

Top 


noting that 
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and that 


II * (n * 7i3) 


30 “ I 30 * 

3r - r 0 ^ -* 


Heace (2.24d) can be decomposed Into the two conditions 


lii + ii 

3z 3r 


G on z 


1 . 


0 r < a 


(3.5c) 


and 


“■(- — -0 on z*l, Osr<a. (3.5d) 

3z r 3i^ 

The lateral boundary conditions (2.25) are 

• 0 3 8w n y / ^ /"I 

— ■ u •• (rv> “ *0 o»n r = a, 0 < z < 1. (>.o) 

3 r 3 r 3r 

In summary, therefore, we are requlrea to solve the equations (2.19)- 
(2.21) subject to the boundary conditions (3.4)^(3.6). With C * 0 and the 
Bond number G no longer appears, Davis and Homsy so we have in 

place of (2.23) the dependence 


M » M (R,Pr,h,a) 
c c ^ ^ ’ 


(3.7) 


for the critical Marangoni nunihor 
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A . LINEAR STABILITY PROBLEM 


The critical Marangonl number at which the conduction aolution loses 
stability Is determined from llnoariz.iClon of the system (2 . 19)-(2 .Zl) 
together with the (llneat) boundary conditions (3*A)-(3.6). Although this 
linear problem Is not aelf-ad joint, we assume the validity of the principle of 

exchange of stabilities, namely that the growth rate of the most dangerous 
disturbance changes from real and negative to real and positive as M increases 

through its critical value. Vidal and Acrivos {,19Gb) show for the linear 
problem on the infinite layer that this is true and so we apply the same 
result to our case. 

When the principle of exchange of stabilities holds, the governing 
equations for the linear stability problem at criticality are 

7^v - Vp + - 0 (^*1) 

7*v - 0 (4.2) 

V‘(j + Mw - 0 (4.3) 

subject to the boundary conditions (3.A)-(3.6). We apply the operator curl 
curl to equation (A.l) and then take the z-compo!ient of t!ie resultant 
equation, we obtain 


..A 

• w ^ 


-I ^ 

M RV^O 


0 


( 4 . 4 ^ 


where Is the planlora Lapl Ian. Thus (4.1) and (4.4) constitute a pair of 
equations for the uaka.;wa ruUv'.lons w .lod r*'e appropriate bounhirv 
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condlClontf are determined by slmpllf/lng (3.4)~^«6); they ere found to be 


3w 

B - w • — -0 on e • 0 , 0 v r a, 

3z 


(4.5) 


+ hO ■ w ■ - 7^0 "0 on z ■ I, 0 < r < a, (4.6) 

Js ... 1 

Oz 


DO 3w . n ^ y y 

■^—■■^■0 on r * e , 0 < z < l . 
or 3r 


(4.7) 


The system (4.3)-(4.4) voduces to a pair of ordinary differential 
equatlora through separation of variables. It was In order to effect this 
reduction that we introduced the artificial condition on the vortlclty at the 
lateral boundary. Ue put 


w(r, ■!), z) - cos mv- (\r)Y(z) 


0(r, 41 , r) - coa m4 J (\r)X(z) 


(4. ft) 


where m • 0, I, 2, ... Is the .tzlauthnl wave number , .I|jj la the Bessel function 
of order m, n.id \ > 0 Is determined by the equation 


J'(Xa) - 0. 


(4.9) 


The condition (4.9) encores that the lateral boundary conditions (4.7) a.*e 
botli satisfied. Substituting (4.8) Into (4.3)-(4.6') we obtain the equations 


- \‘)X t MY • 0, 


(D 


2 2 

X >Y ~ 


-1 

M RX‘; 


(4.10) 
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where D denotes differentiation with respect to z, with boundary conditions 


X - Y - DY - 0 on z - 0 (4.11) 


and 


DX + hX - Y - D^Y + X^X - 0 on z - 1. (4.12) 

Apart from notation, the boundary-value problem (4.10)-(4.12) is 
identical with that solved by Nleld (1964) in his determination of the 
critical Marangonl number for a layer of unbounded horizontal extent* The 
difference between Nield's problem and the one prosencly under consideration 
lies in the significance of the parameter for th.o unbounded layer X is 

wave number In the horizontal plane and can assume all real values, whereas 
for the finite cylinder X is restricted to the set of values defined by (4.9), 

Nleld (1964) solved (4 . l0)-(4 . 12) by expanding In Fourier sine series 
on the Interval (0, 1); complete details of the calculations can be found In 
the cited paper and are accordingly omitted here- It suffices to observe that 
the system (4.10)-(4.12) contains the four parameters M, R, X and h, so that 
nontrivial solutions exist If and only If a functional relationship of the 
form 


MM. R. V. h) - 0 


(4 


1 0 


holds these parameters. Nleld (Ph>4) obtairieii an explicit 

represent at Ion ot <4,13), .aiiioiv 
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M - M(R, h) - W/U 


where W, U are given by 


« d^-R • (nn)^d , . ,v2 . 

W - (1 + h + 2\^ I -^) I 2- + 2RX^ ( I 

n“l n n*l n n*l n 


T “ / ,\2 “ (nTi)'^d T “ / „\2 „ ® (nir) d cosnn 

U . 2X^ ) iHlL-. V n _ (nit) cosnn . j " n 

^ A ^ A ^ A ^ 


, 4 A 

n»l n n“l n 


n*l n n*l 


A » 

n 


with 


2 2 
d i (rnr) + A 
n 


and 


2 

A = d' - RX . 
n n 


It Is convenient to write (^*9) In the form 


J*(s ^ 
m ml 


) - 0 , 


X 

ml 


ml 


/a 


where denotes the i-th positive zero of Thus the integer 

1*1, 2, 3, ... can bo regarded as effectively a radial wave number , 
functional form of (4.14) Is 


M - M(R, h. a, m, 1). 


iiTi Jii !■ 




(4.15) 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


and the 


(4.20) 
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Th€ critical Marangonl number for fixed R* h and a Is then defined to be 

“ M^,(R. h, a) - M(R, h, a, m, 1). (4-21) 

Computations of the Marangonl number as a function of the other 
parameters have been performed using the formulae (4 . U)-(4 . 16) , the Fourier 
series, of course, being suitably truncated. In general It was found that 4- 
diglt accuracy could be achieved with 10 terms. Figure 3 shows the variation 
of M with aspect ratio a for fixed values R ■ 0, h • 0, and for the six wave- 
number pairs m - 0, I, 2, 3, 4, I - I and m - 1, 1-2. These six were chosen 
because for moderate aspect ratios, 0 < a < 2.*) approximately, the critical 
Marangonl number occurs for one or other of them. We see from Figure 3 that 
M - I, I - I Is the critical wave-number pair for small aspect ratios, 
a < 1.15. On the Interval 1.15 < a < 1.65 (approximately) the critical mode 
has m - 2, 1 - I, and with further Increase In aspect ratio this Is replaced 
by the axlsymmetrlc mode m-0, 1 - I on 1.65 \ a < 1.9. Next, the mode 
m - 3, I - 1 la critical In 1.9 < a < 2.3, while for 2.3 < a < 2.5 the modes 
jj - 4^ 1-1 and m - 1, 1-2 give nearly the s.'.me numerical value ol the 
Marantjoni number. 


Figure J 

The ordering of critical modes just described Is retrained for other 
values of R and h. In fact, this ordering Is a direct consequence ot the 
ordering of the numbers Jc! Ined by (4.1'^l, which Is In tut a a consequence 

of tt'.e sidewall boundarv c.ndltlons; It Is nv>t uirprlslng, therefore, that It 
should be Invariant with respect to other ( lyslvul parameters. 
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Figure 4 Illustrates the vsrlstloa of critical Harangonl number 
with aspect ratio a for different value of Rayleigh number R* and at a fixed 
value h - 0. These curves show that decreases as R increases for each 
value of a. Although not illustrated^ computations show that the same 
behavior (M^ decreasing with increasing R) occtirs when h 0. 

Figure 4 

Figure 5 depicts the variation of with a for various values of 
surface Biot number h, and at the fixed value R ■ 0. We see that 
Increases with h at each value of a. Computations show the same tendency at 
non-zero value of Rayleigh number 


Figure 5 


The general pact>?rn of behavior descrlht-d here is consistent wlr’i that 


obtained by Nleld (1964) for an unbounded layer- 
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5 . EIGENFUNCTION EXPANSIONS 


We propose to study the noullnear stability problem by means of a 
modified Galerkin procedure. We represent the field quantities by series of 
functions of the spatial variables, with time-dependent coefficients. 

Following a suggestion of Eckhaus (1965), we shall take as the basis functions 
the eigenfunctions of the linear stability problem. The time-dependent 
coefficients will then effectively be the amplitudes of the appropriate 
convective modes, determined from nonlinear ordinary differential equations to 
which the governing partial differential equations reduce. The series are 
truncated in a rational way, according to criteria discussed by Rosenblat 
(1979). 


Galerkin methods require that the function batils of the expansion 
should consltltute a complete set In an appropriate sense* In this 
regard it can be shown. Appendix A, that Che Marangoni 

nuiuber cannot be used as the eigenvalue parameter on which to construct a 
complete set of elgenf uncc Ions of the linear stability prohlera* This Is 
because the expression gives M as a single-valued fuutlon of R, other 

parameters being held fixed, and corresponding to this there will be only a f^ingle 
eigenfunction* On the otlier hand, if (4.14) Is solved for R, the resulting 
expression has the form 


R - R(M, 1\, a, m, 1) 


(5.1) 


and Is not single-valued; In fact there Is a coviatahly Infinite number of 
solutions of the form (‘>.0 to tt^e equation (4,14). In other wotvls, chore are 
Inflnltelv many values of K fo each value ot M, and correspondi ng ly 


-28- 


infinitely many eigenfunctions. 

Although the critical Marangonl number Is given correctly by the 
results of the previous Section, the implication of the foregoing paragraph Ir, 
that the Rayleigh number is the “true" eigenvalue parameter of the linear 
stability problem. We need to take this into account In setting up the 
Galerkin procedure, even though much of the subsequent analysis will become 
redundant through truncation and approximation. 

Let M be a fixed value of the Marangonl number and consider the linear 
eigenvalue problem 


V^v - Vp + M~Voz - 0 (5.2) 

7*v =« 0 (5.3) 

7^9 + Mw**0 (5.4) 

with boundary conditions (3.4)-(3.6), and with the Rayleigh number R regarded 
as the eigenvalue parameter. Nontrivial solutions of this boundary-value 
problem exist for certain values of denoted where m is the azimuthal 

wave number, 1 is the radial wave number, and J =» 1, 2, 3, ... is the 
particular value implied by (3.1). Thus 


R^ij - h. a, m, 1) 


(5.5) 


and we assume the ordering < ... for other parameters fixed. 

Graphs of the functions (5.5) -an be found In the paper by Rosenblat, lloinsy and 
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Davis (1981), The integer j is in effect a vertical wave number . 

Corresponding to each eigenvalue ^mij there is au eigenvector 

^-mij linear boundary- value problem. The forms of and 

6 are given by (4.8), while the other two velocity components can be 
mi J 

calculated from (5.2) and (5.3). The components, which are required in the 
subsequent computations, are found to be 


u , - (i/\ ,)cos m^ rz) 

mij mi m ml mij 




(5.6) 


w =• cos m>^ J (\ ,r)Y , ,(z) 
mij m mi mil 


3 . ^ - cos m$ J (X r)X (z) 
mij m mi mij 


where X^^j, are the eigensolutions of the boundary-value problem (4.10)- 

(4.12) with R - R_, , , M = M and X = X , . 

J mi 

The explicit forms of the functions X and Y, and their various 
derivatives, are obtained from the Fourier series representations of Nield 
(1964). We omit all details of the calculations but, for the sake of 
completeness, we list in Appendix B the forms of the functions. 

In the modified Galerkin^Eckhaus methi^d to be used below we require 
also the adjoint elgenfunct ions. The system adjoint to (5.2)-(5-4) is easily 
shown to be 


~ Vp* + ^ 


0 


(5.7) 
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7*v^ • 0 


(^.8) 


V^e* + M~^Rw* - 0 


( 5 . 9 ) 


with the adjoint boundary conditions 


0*«qA*v**w*«O on z-0, 0<r<a, 


(5.10) 


+ hO* • w* - “0 on z • 1 , 0 < r < a, 

3z 3z 5z 3z 


(5.11) 


(rv*) * ® 0<z<I. ;5.l?) 

Thesti forms are consistent witii those of Davis (1969) and Davis and 
Homsy (1980). 

The adjoint problem (5.7)-(5.l2) can bo solved by the same Fourier 

series method as the direct problem, and, naturally, the eigenvalues 

the same and have the representation (5.*v). The adjoint eigenvectors are 

denoted (v*. ^ , 6^ ); a relatively simple calculation gives the following 

-mlj *®1 J 

explicit representations: 


u^ - \ cos J’ (a .r)DY* (z) 
mij ml m ml aij 


v^^ ^ • -(o/r)sin J (a r)DY* (z) 
mlj m ml mlj 


C5.Ii) 


:j 1 J ml m ^ ml mlj 


xl ) 


cos ma> J ( . , r )X* (C ) 
m nl mi j 
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wh«re X*j; j, are eigensoluctons of the boundary-value problem 

(D^ - X^)X* + m’Va^v* . 0, (D^ - X^)^Y* - MX* - 0 (5.14) 

with boundary conditions 


X*(0) - Y*(0) - DY*(0) - 0 


(5.r.) 


and 


DX*(1) + hX*(l) + X^DY*(1) - Y*(l) - D^Y*(l) - 0 (5.16) 

for R - ^Ij* forms of the funtions X*, Y* are given In the Appendix B. 

Our purpose Is to study the nonlinear evolution of disturbances as the 
Harangonl number Increases through Its critical value, and at a fixed Rayleigh 
number * To simplify the discussion we take hcncoLorlh 

R - 0 (5.17) 

so that pure Maraagoni convection will be ex.unined. a further 

simplification, and again without essential loss of generality, we take 
the surface Biot number h ^ 0. ITic nonlinear system (2. 19) - (2,21^ can n- w 
conveniently bo written in the form 

' -1 

V^’v - Vp » MPr 't'(v«V)v; (5.18) 

— I? t — — 


7*v - 0 


(5.19) 
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7^0 + Mw - M (|^ + (v*7)0| 


(5.20) 


and the boundary conditions are (3.4)-(3*6) with h ■ 0. 

We shall solve this problem for values of M close to critical by 
expanding the field quantities In series of Che eigenvectors 
time-dependent coefficients. First, however, we let denote an 

eigensolution of the adjoint problem (with h ■ 0) in the case that M ■ the 
critical value at a fixed aspect ratio and with R - 0; the corresponding 
eigenvalue is Because of the ordering of Rayleigh numbers stipulated 

immediately following equation (5.5), this means that • 0 for some m and 

seme i. Let (v, 6) be a vector with div y - 0 satisfying the boundary 
conditions (3.4)-(3.6), and consider the expression 




(5.21) 


M arbitrary, where <•> denotes Integration over the volume 0 < r<a, 0s;< 2n, 

0 < z < 1 occupied by the fluid. Integrating bv parts and noting again that 

(v* , 9^. .) solves the adjoint linear problem with M - and R » R_i we 

-mij mij ^ 

easily find that 


(M - M 




m'^h 

c ml 


0 >. 

ml J 


(5.22) 


Observe that Q ^ 0 when M • and R^j^j " vhlch Is consistent with the 
definition of the linear scablllty problem and Its adjoint at criticality. 

Next, let (v, ’ 3 ) vienote i solution ot the nonlinear system (5.18)- 
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(5.20) with boundary conditions (3.4)-(3.6), for some given value of M. Take 
the scalar product of (5.18) with v*. the product of (5.20) with 0*,., add, 

mlj ffllj 

and Integate over the fluid volume. Using (5.22) we then obtain the equation 


1 36 1 

(M - M )0*. ,w> -MR, ,<w*, ,0> - M<6*, , + Pr v*, 

' c mlj c mlj mlj mij 3t -mlj 3t 


-1 


(5.23) 


+ M< 0 * .(v* 7)0 + Pr V*. .•(v 7 )v>. 
ml J "ml J * " 


Now choose a flaite set ^ of elgensolut ions (v , 0 ) of the linear 

"mlj mlj 

stability problem. Let N > 1 be the number of elements in^, and for 
convenience write 


J - Imljl 


(5.24) 


which means that an element of the set has azimuthal wave number m, radial 
wave number 1 and vertical wave number j. We can thus refer to mlj - p, say, 
as the vector v/ave number of an element 


We assume that the solution vector (v, 6) can be represented, to a 
good approximation, by a linear combination of elements of the sec^, with 
time-dependent coefficients. Thus we set 


(v. 0> 




(5.25) 


substitution of (5.25 into 0*23) reduces the latter to a system n N 
ordinary nonlinear differential equations for the amplitude functions 

The details of tn! reduction are cons lder.jbly simpliflri I'n a:* )ani 
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cf the following orthogonality relations. First cf all we tU<-? 

bi-orthogonality condition 


<w*6 > » 0 when p i* q (5.26) 

P q 

for any two fields with vector wave numbers p, q. Next, since the azimuthal 
dependence has the form of trigonometrical functions, we deduce from the 
latter's ortlio^'uallty properties Chat 


''^*4 i'' I 0 ^ “ I when m. ^ n (5.27) 

oij nkl -mij -nki mlj nki ' ' 


foi any values of i, J, k, i. Similarly, since the radial dependence has Che 
form of Bessel functions, we have that 

, pivGt 

. ,-p< <0* ,w , > • <v*^ <6*. ,9 , a> ■ 0 when i k (5.28) 

ri’H'-'* mlj nk£ -mij -nkl mlj nki 


for any values of m, j, n, t. :T- 

There remains the question cf the choice of the set^ . For ease of 
computation it is desirable that^^ should comprise as few elements as 
reasonably possible. Next, since we are concerned with the weakly nonlinear 
interactions that cause the onset of convection, we must certainly I u lui'.e I n ;<{ 
Che critical mode (or n ) .is pointed out In ihe or-. ; >as Section the 

nature -if i h.- i rftlt ‘1 mode depends on the aspect ratio. Finally must 
include a minimal number of other modes ne. I ’ i . ^jinerate no.nlinear evolution 
Into convection of the critical mode, hy “niin nAr' we mean the non-crltlcal 
modes t>i th.- '-•i l’i.gt damping rates. As can be seen from (5.23), when M Is 
i J the damping rate of a non-crltlcal node Is determine’ 


ipl.roxlmacely by the magnlta.|.> .r so that we retain only those -nolos 

wtih the smallest values of j and neglect all oth,.,-,, \ i; ;.-,,s8lon of the 

...lore and Its validity can be found In Ro:.enblat (1979). 
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>.VO<.i{TT0N AT SIMPLE P OINTS 

We study lx ' S.-. Uou weakly noallnoar evolution Into convectlca at 

three specific valva*-i •»* i’m’ aspect ratio, namely a - 0*90, a *l.SO ai\i( .. * 
l,510* As can be seen from Fl^jure I, the loss of stability of the basL^ 

* • utnctlon state at each of these values Is s i In the sense that only one 

mode loses stability as the Marangonl nuuh^ I . through its critical 

value. In each casi- i» the problem to a single ordinary nonllueic 

d ( f feront Ivil equation of Landau type, and examine its solutlon^i ‘ id i r 
st dll llty . 


A. The case a • 0.90. 


Figure 3 shows that Che crlttcal it ti)ls aspect ratio is the mode 

in, that Is, azimuthal, radial and verticil \umhcrs all equal to 

uMtry. We find that 


- 79.5 (b.n 

and, by h/oMi lies I s , ^ - 0. Now the quadr if tc sol f-lnteract ion of the mode 

III generates the modes 0(J, wii)\ ) * I, 2, ... . Computations show, 

however, that 


"'-HI 


l.J 'oij* 


nln ^ 

^211 ' 1,J ^21J^ 


^0.2) 


and, morc 'Vi'!', rh i* ; ’i.»sc Rayleigh numbers are widely separ.ii» ‘ '\[i^ 

reason we t iL ‘ ;* - t ^ ’'iprlse three modes: 




j’ - {ill, oil, 


(6.3) 


and substitute Into (5.23) t'n* 


(V, 0) - 




^11> 


^011^-011’ 


Vi> 


'^211^-211’ 


®2ll>* 


(6.4) 


where the 


are functions of time. 


Using the orthogonality relations (5.20)-(5.28) to eliminate several 
of cU.* lemiSjWe obtain the following set of equations: 


''lll^Ul * ■ \^^111 ' ^111 


- 1 . 


''OII^OII ” (M - '^Oll^OU^^OU ^011 


^211*^211 “ ^’c " '''c ^211^211^^211 ■ ^211 


(6.5) 

(6.6) 

(6.7) 


where the prime denotes differentiation with res. to t, whore 

(6.S) 
(6.9) 

11 .: w'l-r.* the ari hotnogoneous quadratic functions of and 

The genera) forn of r'l' 2^,^ is g’von by 


'i'0*. .9 . ct-?r“^v*, *v . ,> 
gj j ml j j j 


mi J 


4« 4 4> 

mij mlj 


jn' 


<W*4 4^ 4 4> 

uil j ml j 

^alj " <0*, .w , ,> 
mij I'llj 
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<8*v >Z 
P P 


P 


M<9*(); A 


V 

q q 



A 0 
r r 


♦ Pr~^v*< 
-P 


(i A V •7)y A V > 


where p, n, r r<'f. r to vector wave tuunbers. 


( 6 . 10 ) 


Each of the two terms on the right >' ('..Ih) contains Integrals 
when^ has N elements. Mi-w.-v r, -levcral of those Integrals are Identically 
tero by virtue of the orthogonality rel>c^>n^ (5.22)- (5 . 24)* In the present 
case we find that, f<>i- .’i.> .ooje 111, only 4 of the 9 Integrals 
are non-zero, and we obtain for the as.., t .i.l .[.mlratlc nonlinearity ai 
expression .if i'n« form 


III Ul 


^'lll^ V'^Oll 


^ 2 '' 21 I^ 


( 6 . 11 ) 


where 


d 


mlj 


<e*, W > 

ra I J i»i I j 


(6.12) 


.ui(l whi re 




lll'--lir''’all -tall 




'I 


•( V 


Vv 


- Ill "-lil -la: 


f V 


-T1 


►Vv )> 

- - 1 1 r 


(b.l 


for a • 0 .ind « ’ '>1 'iK^r .'onsldvn .i? ions le obtain 


Oil 




ii 


■'oo' 


i > ' * 


'..l' 


(t.U) 
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‘* 211^211 




(6.15) 


where 


mk 


< 3 * / 

mll-kll 


VO 


kll 


/ • C ' ' 

-rail ^-kll 


•V)v > 
''-kll 


(6.16) 


for m “ 0, k ^ 0, 1, 2 .uiil ra ■» 2, k =» I, and 


“202 “ "^*^ 211^-011 *'^^211 '*■ - 211 *^^ 011 ^ ^ - 211 *^-' 


• Vv 

oil -21 




(6.17) 


Equations (^» • • 7) are the evolution equation for the mode 111 at 

the aspect ratio a • 0.9. The null soluiton, » ^211 * 

oorr.--;iouJs to the conduction state, and at fixed RayloI^N nuiahor,R ■ 

Is stable for M < and unstable for H > 


To study bifurcation from the critical the evolution of 

iouv.-v :lt>n in the neighborhood of M • M^, we can simplify sy.i^ i (').5)- 
(6.7) in the following wiy. T'\ ‘ nodoi; .\q|| and A 211 relatively strongly 

damped at M » and are pr.-..* U i.nly due to the quadratic self-interaction 

of the >*ri!i\.il node Hence , we can neglect the t ime-der 1/ it Iv.* i.-riw i 

(6.6) and (6.7) and replocf M hy f »i these equations. Moreover, when M Is 
close to i'i»‘ iiignltudes 01 ^Qii ^211 snuill o>ap.4r«*d with Ai^|; 

Hence in the right-hand >1 >r (T-.H) and (6.15) we can neglect ih.- 

quadratic terns involving - \‘l A^^i hy c^)mparl8on with the terms l ic.lvlag 

^lll‘ Takiny, Ihfse approximations together, and substituilug (6.1.) i>id 
(6.15) iu> (h.6), (6.7' respectively, we obtain 


dll 




M*“a, 


; I 


\\\ 


Kid 
211 2Ll 2U 


(6.13: 


Mi 
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We now substitute (6.18) Into (''.n)! > < ^ iheti substitute the latter, 
with M replucfil hy M , !, ito (6.5). This glvos th.- single equation 




(6.19) 


where 


(j 


111 


M 


■ill 




^2 '*21 


^01 1^01 1*^01 1 


’^211^211*^211 


--•). 


( 6 . 20 ) 


Equation (o.r?) I- t’i<- T,.<ndau equation for the evolu' Ion o*"' tlic 
critical mode III. The coef f i . I »*m ^ '^111* '* ’etermined by numerical 

Integratlo I .f tli.- t.>i>iopr late prodici ■. >r ’ifunctlons. The computat Ion-; 

have been performed at various . of Prandtl number and some . i>’.- 

-.ii.'Oi in Table 1. We Infer from the calculations thi' 


111 


> c. 


Ill 


> ■) 




0 ! 


'1 


fjr »I ' "i ri ltl 'uimbers. 


Vr 

‘111 ^ 

"ill • 

O.l 

0 . T / 

1 • •- 

1 .0 

O.IJ 

0.. lO 

.0.0 

O.IO 

O.ll 


O.IO 

0, 10 


O'- : Wi. 1 ;,G£ IS 
OF FCCH QjALITY 


Table 1 


-41- 


Prom (6.19) and (6.21) we infer tli.ii < blfuractes f rimi Llic* 

critical value having the repreaantti l<>u 

^11 * ^ 

The aolutlixi only iior M>M^ (supercrltlcally) aad knowii from 

elementary Mfurit itlon theory to be stabL ^ Because of the rnpr.-*.. o itlon 
(6*4) |We conclude that at aspect ratl*i < ih.- onset of convection is 

supercritical and^ r..r ^ lall M - has to leading order i h«- r.r.i ^ non- 
.i\f sytiimetrlc raode with azimuthal wave number 1. 

The case a - 1.50. 


From Figure 3 we see chai at (hi . ^apect ratio the critical nod-> fs 
the azimuthal wave number is 2^ tud the ridial and vertical wave naohers 
are both unity. The critical MiraJig-»n( aiahar Is 


“ 7^.5 (6.23) 

with ^ “ 0. The quadratic Interactlo i if - ii . ’ . with itself 
genera':.*- ' )i.> ), 4iJ, with 1, J - I, 2, T '..i mse of the 

irderlng .>f 'J-.e >. i iiej Rayleigh numhcrs we approxl aot •• ’>. .•,,l.‘:tlng all 

except f'. -n 01 ; inrl -til. Thus ^ 


-42- 


J - {211, OU, 4!lj 


(6.24) 


aiul 


(v, 0) “ ‘^21l''-2ir *^211^ *’ ‘'oil^^ir ^Oll^ \u^-4ll’ 

We substitute » ii»> (S.23) to v)btaln a system of three )^*t>aary 

differential equatloa-- r,.r i ’le amplitudes, namely 


''’U'2U ' ~ ^^''211 ~ ''211 


(6.26) 


^ 11^011 " ~ ’^Oll^Oll^'^Oll ~ ^011 


(6.27) 


v,,,A,,, - (M - M - M~^R,,,f,,,)A;,. - 2,,, 
411 411 c c 411 411 411 4il 


(6.28) 


where the coefficients a ic * iht‘ forms (6.8)“(6.10). Proc.v !l as In the 
previous w<- f 1 id that 


'^211^211 * ^2ll^^o\ll ^4^411^ 


(6.29) 


where 




”.ll ^ S.11 


■’® 2 U> * 


-211 




V V ^ V • r ^ ^ 
211 -mil ^11 - 2 ir 


(6.30) 


f <r ii •* 0 .iiul :ii ” ; ilao. 
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(II 


^011 “ ”^*^ 02^211 *^ 00^011 


^ 3o4^11> 


(6.31) 


nriK 


‘^411^411 


”^^42^211 '^404^01l\u^ 


(6.32) 


where 


-1 


2 ;n-< 2 vi,-''>Svii> ( •”> 


mil '-kll '-kU 


for m ■ 0, k 0, 2, 4 <iii.I m ■ 4, k ■ 2, and 
■ <“Ju'i0U-”4U " -4ir%U> * 

(6.34) 

Using the same reai>iU i;:, i'- hvfore^we can solve (6.27) and (6.28) 
approximately to ftn.l ‘v). . lod ^ 41 ^^ In terms of Aoj]- Wo olo ,ln 


-M ^ ^ 

‘ ■ “ ,C, , ■ 


*^42 ^2 11 


(6.33) 


-.1 rm 411 


Subsc ! I It i (6.35) and (6.29) into (6.26), we f^i* ’i.tt l ht- latter reduces to 
the simple Larului e r i • i I o i 


^ 2 ir ' 21 l ■ -^^^211 ~ '* 21 l '^211 


( 0 . 36 ) 


where 


“^o ‘^02 


-♦ ‘4 2 


> 1 ' 


'2 ' '^n. *'i ."HI '^^ll^41l‘^411 


). 
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Compaf:<' 1 v J i.--s .>f ttie coef flcl'ini > ^211’ ^211 Prandtl 

ir«* jjlven In Table 2, from whl :h li > • ’>• ^.'t-n Llut both coefficients 

are always positive- We Inri-r '’it' 

ri-prostints a stable supercritical coiiil i i i * i ^ i^nlion for M - small, and 
. u r,-sponds to a aon-axls>« » u i r* J . V)de with azimuthal wave number 2. 


-4 

Pr X 10 X 10 


0.1 


0.37 6.2 


1.0 


0.13 0.98 


10.0 


0.10 0.50 


o> 


0.10 0.45 




Tabic 2 

:"GE fS 




‘ > C . ALirv 

C. The case a « K80. 




Figure 3 ; ii 

‘ he 

critical mode at this aspe:t t <M 

1 . . ' . ‘ 

axl ‘'•’'Tomet ric mode Oil. We 

f ind 

th i' 
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with RqiI " quadra: !>• r-i i»i e-ractloa of this mode generates ill the 

mden OlJ with 1 J ■ 1, 2, but by vlrim- •T i*. irlerlng of the Rayleljjh 

numbers Rqij retain only the mode O'l’ . T'l i- 




- (oil, 0?1 




-tiid w<* suhetltute 




Into (5.23) to obtalri i ,>it( if ^iinplltude equatloa^i: 


'^01 ” 'o^^oii ^011 


(6.42) 


^)2l'\')2l ~ *■’’ ~ ’c " \ V21^02l^^021 “ ^021' (<’•^3) 


i'sin« ‘h- f-H-Qula (6.10) we ftiul (.li.ii 


*^OU^iUl “ Oil ^U2''oil''o21 'l2‘'02l^* '• ' 


'■^ 021^021 " ^^^ 21^011 ^ 212 ^ 0 ll'^ 02 l ^ 22 ^ 021 ' 


r M . 4 •) ) 


. rn 


Y =■ ' ★ f ^ 9 P r * • f V • ^ V ^ 

’iK Oil' -Ok I ■' Ok! - II : '' -v..: 


'■>) 


for I, k ■ I, 2, .in' 
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-1 


’lU ■ <®3u<Soii'''V" ^ ion ''2011 -’iozi “ io2r’ion>^' 

(ft. 47) 


We approximate as befori- rut* M < to and on the assumption th.it 

th«^ f ade of is much smaller than 'r We then solve (6.43) 

and (6.45) to ahtai.n 


^021 “ R 


-''^ 2 l*on 


02 1^02 1 *^02 1 


(6.48) 


Wh substitute this and (6.44) into (6.4?), vlih M rt-ilaced by in the 

2 

uotilf \.- ir terms and with the term c.)mI ilulng Aq 2 i omitted on the groan 1^ Ihni 
l! Is smaller than those retained. This 1 to a single equaitons *:or i 

i r f t leal-mode ampll tu<le, y 








49) 


where 


Y 


0 




'oil 


(6.50) 


and 


"Oil 


>w *. i i 

lift,, d ■ . 

) 1 ; ill! 1 1 > ' II . 


( 6 . 51 ) 


0 .•nput 'd valj'ii ’ • [i-nts are given In 3. 
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Pr 

1 

0 

X 

1 

0 

> 

Y„ X 10-^ 

uj X 10*^ 

“oil 

0.1 

0.36 

0.32 

1.8 

1.0 

0.12 

-0.19 

0.23 

10.0 

0.10 

-0.24 

0.15 

go 

0.98 

-0.25 

0.14 


*V\h 1 ^ 

Observe Chat hoib .v.islMve, hut Chat I 9 positive foe 

I.w ''i i.ulil numbers and negative for raoderar:n ^ Prandtl numberj- 

One solution of (6.49) Ls ^ - 1 ), vhich corresponds to the 
cuiibniion state. This solution Is stable for '{ ' M ii«! nusL.ihle for 
M > Ocher ir ’ lucBlned from roots of the equation 


‘^()u Sni '^o 'oii 


- t 




0. 


(^.S2) 


AU t>i“ 'Solutions are Illustrated in Figure 6a ^ 0, and In 

6b f'^r the case Y.^ - A conduction solution t^..l .. - for both M < 

1 vl M ^ •'Lclcal bifurcation), but t’\e subcrltlc*! br.in h ^'ns 

-trmind at a value Mq of M and contl ii - ( i' * fu* half-plane M > 

'‘^oeis 

9u,Aurv 

■• ’•'.•I- il l. I'lo -it-ibilUY Jf the sol.itl>\ ' ■ '.’Inwini^ way. Let 

ienCe any tl’i -i ’ dilution of {k . ~e set 


OWGW/,1 

O'^POOfi 


In equation (6>49) and linearize, to obti<-i iln- hI ability equation 


*0 ■ - ••‘c • ^’ 3 * • Nu*'>‘o- 

A 8l«;>l-^ t1-Ml itlon -ibiws t'lai Mt." aiiWrl t leal branch OP In Figure 6 Is 
unstable, and tl*.at br.neheo OR, PQ are stable* This Is the 
result '•»* t ranacrltlcal bifurcation* 

These calculations settle In j»rlu.-!,il • .tie question of the direction 

of the flow at the cuuf ~f ><r 1 %•* «-outalner because .tsynnn«try of the 

bifurcation diagram Impl t *s a ;ireferred branch* As M lncreaH«*s l•»v.l•*is M^., a 

disturbance, <mall, to the conduction sol ii i<>u wfll result In loss of 

stability of the lati--r‘ •x-Poce M reaches lli«- «rl ileal value M^, and a 

consequent snap-thr ough to Ch-* branch PQ. (See, for exaaplu, Rns-nhlat 1979. 

for a riin.-r .It >«i*.a 38 lon of this process*) As M Incre.ts.-. ^t 111 futher the 

system stay', "n the branch PQ which Is thus "preferred" !•» Ihe branch OR* 

Of course, when M > M one could find a disturbance of large enough size 

c 

that would cause the system to Jump from this branch to the other. The 
fact that both supercritical branches are stable is similar to the result 
of Liang, Vidal and Acrivos (1969) for axisymmetric buoyancy driven convection 
in a cylinder. 


Ch J.U. h/iGE <S 
OF k>00H gUALiTY 



From (5»(>) Wrt 


c'ut th« vertical vtloclty coaponcnt f»r ««t*t 


Oil tj 


t- 


Oil 




(6.5S) 


-ii i'i«- of Ch« contain****. ‘JfwrlrMl .'.tlculaClono show that ^oil^*^ ^ ^ 

0 < a < 1; hsrs tha sign of Aqj^j d«t**rnl .;.*•* i*i*» dlractlon of the fl<»w, 
giving upflow whs.i ‘Vh ^ ^ ^ '*'0 'lownflow %rtton K 0» Tha prafar**»»l 

PQ haa < 0 f.»r vary aaall Prandtl nunbcrn anJ > 0 f.n* wiittraca and 

Prandtl nuabars. He Infar that that** vM! it.* .hiwuPlow at tha cantar 
w’inn Pr « 1, which la tha cav«> for ll*|ut t >i**tilst upflow whan Pr > 1, 
which for coianon liquids. 
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7. EVOLUTION AT DOUBLE POINTS. 


Ic Is evident from Figures 3-5 th.tC there are certain values of Che 
aspect ratio at which two modes lose stability simultaneously. In this 
Section we shall Investigate the onset of convection In Che neighborhoods of 
such double points and the secondary bifurcations that can result. 

A. Intersection of modes 111, 211. 

Ac the point marked A In Figure 3 Che curves of N as a function of 
aspect ratio for the modes 111 and 211 Intersect. The value of a at which 
this Intersection cakes place la denoted a^^ 1.20). The comnon 

value of at this point will be denoted computations jlve 

”c " ”iii “ ”:u * 

By hypothesis we have that Rm ■ ^211 " 

We are Interested in studying the onset of convection at values of a 
slightly less and slightly greater than a^^. As noted in the previous 
Section, the self-interaction of the mode III generates modes with azimuthal 
wave numbers 0 and 2, while the self- interact ion of the mode 211 generates 
modes with m - 0 and m ■ 4. In addition, the interaction of modes 111 and 
211 generates a mode with m • 3. For reasons indicated earlier, we approxi- 
mate by retaining or'.y the leading (in the sense of Rayleigh-number ordering) 
member of each of these generated sets. Ihis leads us Co select a S-elcment 
sot of expansion functions, namely 
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(111,211,011,311,411} 


(7.2) 


Substituting the appropriate eigenfunction expansion into (5.23)^ tre 
ob .ain a system of five ordinary differential equations for the amplitudes. 
These equations are conveniently written as follows: 




-^.11 • 


m ■ 1 and m • 2 


(7.3) 


and 


''mll^nll 


(M-M~^a ,,f ,,)A ,, 
c roll roll roll 


roll 


m ■ 0,3 and 4 


(7.4) 


Note that in (7.3) we have retained in place of the reason for 

this will soon become apparent. 

We reduce the system (7.3)- (7.4) as in the previous Section. We use 
(6.10) to calculate the quadratic nonl inearl ties in (7.4) and find 

‘^011^011 * ^^'^01^11 ‘^02'^11 * '^00^11 ^3^311 ^04\ll^ (7.5a) 


‘^311^311 " ^Si2^11'4i1 


(7.5b) 


and 


‘^411^411 “ ”^\2^ll ^404^)ll\ll ^413^11^11^ ’ 


(7.5c) 


where the are given by (6.16), is given by (6.34), and wliere 

°mjk “ ^®mll%ll*^kll ^11*^®J11^ ^il'^jll*^^!! 




(7.6) 


We now solve (7.4) approximately by neglecting the time-derivative terms, 
replacing M by M^, and neglecting as relatively small the quadratic terms 
involving 4^]^^ ‘‘nd on the right-hand sides of equations (7.5). 
This process gives the approximations 
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■ »ou'on“ou 


, •*'^«2'^u 

^4ll^4ll‘*411 


(7.7) 


. , **C^312'S.1I^I1 

^Jll *^3ii^3uSll 


These fornmlae can be compared with (6.18) and (6.35), it bel’«g noted in 
particular that a ^2 ■ Pqo 0 La 9 “ ^a*>* 


%2 '^42* 


Tlte quadratic nonlinear! ties in (7.3) are found to be given by 


‘* 111^111 “ ”^11 ^ 123 ^ 211^11 


(7.8) 


‘^211^211 ” ”^‘^21^11 ^0^11^)ll ^4"^ll\ll t ^l3^11^Jll^ 


where a^, are defined by (6.13), ou,j^ by (6.16), 3^, by (6.30), QL 213 » 
0^23 (7.6), We now substitute (7.7)-(7.9) into (7.3) to obtain the 


following pair of equations: 

''ll All * ^”*”lll^\ll ■ ®l\ll^ll ■ '^1^11 ■ '^iSll^ll 


(7.10) 


''211^11 “ ^”"^11^^>11 " ‘^2'Hll ■ ‘^2^11^:11 ’ '*’2^211 (7.11) 


khere the coefficients are defined by the following formulas. 


• \^2/dlll ’ ' \^l^‘‘211 


_£_ s ~"0 02 ^123' 312 I, 

111 I ^0ll^0ll‘*011 *^311^311*^311 ' 


’o^)l 


“211 [ *^011 ^01 1*^01 1 '*3ll‘31l‘'3ll j 

*^01l^01l‘*01l‘*lll ’ 


^213^312 
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where «>21I (leCineU by (6.37). ConpuCed values of Che coefficients are • 
given in Table 4. 


Pr 

v^xlO-* 

c^ X 10'“ 

X 10’^ 

o^ X 10*^ 

''211 ^ 

C2 XlO’^ 

X 10“^ 

u>2 X 10 

0.1 

0.49 

-1.2 

14.0 

2.0 

3.7 

1.3 

1.6 

4.2 

1.0 

0.19 

1 

o 

• 

1.5 

0.27 

1.2 

0.36 

0.28 

0.62 

10.0 

0.16 

-0.69 

0.93 

0.18 

0.91 

0.26 

0.20 

0.31 

m 

0.16 

-0.69 

0.88 

0.7 

0.89 

0.25 

0.19 

0.28 


Table 4 

Ue propose to study the nature and stability of solutions of (7.10)- 
(7.11) In the neighborhood of the double point A, and for values of M 
reasonably close to M^. 

We consider first the case of aspect ratios slightly less than a^, 

d fC a,. Vie then have that 
A 


”U1 ^ ^211 


(7.13) 


and we define 




(7.14) 


Rquaclons (7.10)-(7.11) can now be written In the form 




(7.15) 
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where for the sake of brevity the notation has been simplified In an obvious 
way. 


Because of the approximations used In their derivation, equations 
(7.15) can be regarded as valid only for small values of A and small values 
of n. A reasonable measure of smallness Is the .'atlo of these quantities to 
the critical value of M, which Is about 85. In the computations to be 
described presently we have taken values of A In the range 0 < A < 5.0 and 
of n In the range 0 s n < 2A. Results for significantly larger values are 
not easily substantiated on ttie basis of our approximation scheme. 

Before proceeding further we formulate the stability problem 
associated with the system (7.15). Let (A^, A^) denote a solution of (7.15) 
and set 


> A^ + a^, A2 - A2 + 82. 


(7.16) 


Linearization of the equations with respect to the disturbances leads to the 
linear system 

• _ _2 -2 - 
'^1*1 ■ " *^1^2 ~ ^*^1^1 ~ ®l^^“l ~ ^^1 

(7.17) 


V2^2 " “2(c 2 + 02A2)A^a^ + (n - A - o,A^ - 3 w 2 A 2 )a 2 * 


The exponents of this system determine the stability of the 
solution (A^, A2). 


We exa*nlne now the solutions of (7.15) and their stability. Observe 
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flrst that (7.15) has the trivial solution ■ A 2 * 0 , which corresponds to 

the conduction state, and which is stable for n < 0 and unstable for n > 0 . 
Moreover, for this solution one stability exponent changes sign at n * 0 and 
the other at n ■ A, so that each of these values locates a bifurcation point 
for the appearance of a new (convection) solution. 


The trivial solution is unique for r> < 0. As n increases through zero, a 
pair of nontrivial solutions emerges, determined by the pair of equations 


2 ■’ 
Ai - (n - Cj^Aj - 


(7.18) 


“2*2 * 




0 . 


If L were large, these equations would reduce to A 2 * 0 and (equivalent to 
^lll^ given by ( 6 . 22 ), and would correspond to a **pure** convection state with 
azinuchal wave nunber m • I at leading order. When a is small, however, the 
mode m ■ 2 has an effect, as expressed by the coupled equations (7.18), in two 
ways: the solution has 9 ^ 0 and the parabolas ( 6 . 22 ) ore distorted 

as n approaches A. 

Numerical calculations show that solutions of (7.ld) exist only on an 
Interval 0 < n < n^, where > A. The value of Oj, depends on Prandtl number, 
but computations reveal that 1 n,j/A <1.2 for the entire Prandtl number 
range. On 0 < n < there Is precisely one root A^ of the ctiblc in (7.18), 
and a corresponding pair of roots These solutions are Illustrated 

in Figure 7: one pair of solutions is represented by the curves O^UTj^, ^2^2 

and the other by the curves O^LTp 02 T 2 « 
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Flgure 7 


Calculations utilizing equations (7.17) show that these solutions are 

stable. 

The remaining solutions of (7.15) are given by 

Aj = 0, n - A • *^2^2* (7*19) 

These exist for n > A and are the "pure*' m • 2 mode solutions considered In B. 
of Section 6. Their stability Is determined by substituting (7.19) Into 
(7.17), which gives 

'^1*1 ■ ■ *^1^2 ~ ‘^1^^*!* ''2*1 ’ ■ A)a2- (7.20) 

Evidently It Is the first of these equations that decides stability. 
Calculations show that the branch corresponding to positive A 2 , the upper 
branch AR In Figure 7, Is unstable, while the lower branch, AS, Is unstable 
Initially but regains stability at precisely the value n ■ defined above. 

The situation depicted In Figure 7 has been established numerically. 
In summary, the behavior of the system Is as follows. For n < 0 (N < 
the conduction solution Is stable, and is replaced on 0 < n < flp (^111 ^ ^ ^ 
My, say) by a pair of mixed-mode solutions given by (7.18). For T] > T|^, 
however, these give way to a single solution which at leading 
order is a pure convective state with m - 2. This behavior is quslitatlvely 
independent of Prandtl number. 


For aspect ratios slightly greater than a > a » 

A 


we have 
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^211 < ^ 111 * 


(7.21) 


Ue then set 


M - M 


211 ’ 


A - 


(7.22) 


to obtain the system 


• 3 2 

* (n * A)A^ * c^A^A2 * **^1^1 "* 


• 2 2 3 

^ 2^2 * nA 2 C2A^ * 02 A ^ A 2 ” **^^2 * 


(7.23) 


The stability equations for a solution (A^, A^) of this system are similar to 
(7.17), except that the terms n and n " A are Interchanged. 

Figure 8 


The results of numerical calculations for (7.23), Illustrated In 
Figure 8, are as follows. The conduction solution A]^ ■ A 2 " 0 exists for 
all n, and Is stable for n < 0 and unstable for n > 0. The solutions defined 
by 


- 0, - Aj ■ n/-*2 (7.24) 

exist for n > 0. The lower branch O 2 S Is always stable, while the upper 
branch Is stable for 0 < n < n^, where c (0, A) and has a vilue dependent 
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on Prandtl number. The solutions determined by 


Aj - (n - A - 

<** 2^2 ^ ^ 


(7.25) 


exist on the interval < n < A. There are two such solutions, comprising a 
single root A 2 and corresponding ±A^. These solutions are found to be stable. 

Ue see from Figure 8 that the conduction solution is replaced by the 
pair of solutions (7.24) on 0 < n < n.^. Next, there is a region, n,^ < n < A, 
in which there are three stable solutiona: the lower branch of (7.24) and the 

mixed-mode solutions of (7.25). However, for n > A there is only one (stable) 
solution, which is a pure nonaxisymmutrlc mode with m - 2. These results 
apply at all Prandtl numbers. 

We see from Figures 7 and 8, and from the preceding discussions, that 
on either side of the aspect ratio a^ the system eventually attains the same 
state: a single convection solution with m > 2. In one way or another the 

mode with m * I is suppressed by the interaction when A is small. This 
interaction can therefore be regarded as being a mechanism for wave-number 
selection in the sense just described. 


B. Intersection of modes 211, Oil. 

The point B in Figure 3 is the intersection of the curves for the 
modes 211 and Oil. We denote the corresponding value of a by a^ (^ 1.70). 
The coosDon value of ^211* ^011 this point will be denoted by 



computations give 
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”c “ **211 “ **011 “ 

By bypochesis ve have chat - 0. Ualng Che same reasoning as 

in Che previous case, we cake 

i • l2U, Oil, All, 021} (7.27) 

as the sec of eigenfunctions for Che evolution of Che modes 211 and Oil near 
Che double point. 

Substituting Che associated eigenfunction expansion Into (^>23) we 
obtain, by analogy with (7.3) and (7.4), the following ordinary differencial 
equations. 


''mll^mll 


- (M - 


**mll^\ll 


- Z 


mil* 


2 and m ■ 0, 


(7.28) 


and 


V ,,A , - (M - M - m”*R ,,f ,,)A ,, - Z ,, 

mil mil c c mil mil mil mil 


(7.29) 


for m ■ 4, 1 ■ 1 and m “ 0, 1 ■ 2. For the quadratic nonllnearltles In (7.29) 
we use (6.10) to find chat 


®4ll^4ll ■“ **^^42^211 ®404^0ll^4ll ®40 a\)?i\i1^ (7.30) 

where given by (6.33), by (6.34) and ^>7 (6.34) with the mode 

Oil replaced by Che mode 021, and 


®02l^02l " **^^22^211 ^2l‘'oil *■ ^Zl'^Oll * ^212^011^211 ^22^021 ^24^411^* 


(7.U) 
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wh«r« Che ere given by (6.46)-’(6.A7), and where we define 


-I 




for k ■ 1, 2 and m 2, 4. Ue now aolve (7.29)t making Che same 
approximacions as before, Co obcain 


-M 


'021 


*^ 021 ^ 021^*021 


2 2 
^*22^211 ^21^)ll^’ 


'411 


-**c^42^fl 

*411^41l‘*411 


(7.33) 


The quadracic nonlinearicies in (7.27) are found Co be given by 


‘* 211^211 


**^ 211^^011 


^ ‘"0''021 


‘^4^11^’ 


(7.34) 


where 0^ are given by (6 >30), and 0^ is also defined by (6.30) wich 011 
replaced by 021, and 


^011^011 “ **^^11^011 ^112^011^021 ^12^021 ^12^11 ^14^11^ (7.35) 


where Che y’o are given by (6.46)-(6.47) and Che 6's by (7.31). 

Ue now subsCituCe (7.33)-(7.35) inco (7.28) Co obcain Che following 
pair of equacions: 


• 2 3 

''211^211 “ " ”211^^211 ■ ‘^2^211^011 ■ V2ir'oii " “2^211 


(7.36) 


• 2 2 2 3 

''OII^OII • " ^OU^^Oll ■ ^O^Oll ■ *^0^211 ’ °0^0ll^211 " ““o^Oll 


where 
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■"r^o’'2i 


021 021 021 2U 
.3 - 


"O R 


‘‘‘c ^22 ^11 2 


021 021 021 Oil 


. ^ ^ 0^22 ^ 

**2ll ’^02l‘*02l^02l *^411*^411^411 


(7.37) 


where defined by (6.51) and where l3 gWen by (6.50). Computed 

valueo of Che coefficients are given In Table *5. 


Pr 

''211 ^ 

c^ X 10‘^ 

-2 

o^XlO ^ 

(4)2 X 10‘^ 


Y„xl0-^ 

'0 

u)qX10 

0.1 

0.40 

0.22 

6.9 

11.0 

3.6 

0.36 

-19.0 

1.6 

1.7 

1.0 

0.15 

4.3 

0.55 

1.5 

1.2 

-0.21 

- 1.4 

0.22 

0.21 

10. 0 

0.12 

4.7 

0.30 

0.98 

0.96 

-0.26 

0.40 

0.15 

0.14 

OD 

0.12 

4.8 

0.28 

0.94 

0.94 

-0.27 

0.57 

0.14 

0.13 


Table 5 


We proceed with the analysis as In the previous case. For aspect 


•atl03 sllithtly less than a-, a < a_, we dcflt 


-' ■ "on - ":u > “• 


(7. .IS) 


whereup equations (7.36) becuire 




'* 2^2 " ^2 ” ® 2 ^ 2^0 ~ ® 2 ^ 2^0 ” “ 2^2 


(7.39) 




(n - a)Aq - YjaJ - 


'=0*i • 




“ 0 ^ 


with an obvious abbreviated notation. The stability problem for a solution 
0 ^ 2 * ^0^ (7.39) Is determined from the equations 

'^ 2*2 " “ *^ 2^0 ” ® 2^0 ” ^“ 2^2 ^*2 " ^*^ 2^2 ^® 2 ^ 2 ^ 0^*0 

(7.40) 

# — — • — —2 
''0*0 " ”^^*^0^2 ^ ^®0^^2^*2 ■*■ (*1 ■ " ^7qAq - 3wqAq)Sq. 

Because of the proximity and flatness of the curves of Figure 3 in the 
neighborhood of the point B, we estimate that the range 0 < A < l.O is 
reasonable for A, with n not larger than 2A. 


Figure 9 


The solutions of (7.39) and their stability behavior are illustrated 
In llgure 9 for Pr ■ 1.0, 10.0 and •; for the case Pr • 0.1 the graphs of Aq 
against n have to be reflected with respected to the n**axls. The conJuctlon 
solution A 2 ■ Aq • 0 loses stability at n ■ 0, and the solution which 
bifurcates from this point Is primarily an m * 2 mode, modified by the 
presence of a small m * 0 component. The latter Is due to the fact that 
(7.39) have nu nontrivial solutions with Aq ■ 0, owing to the presence of the 
term cqA^. However, the modification Is small because cq is a small 
quantity. This solution Is represented by the curves OM, ON, In Figure 9, and 
derives from the equations 




■* 


MMMBitMtUSlimA 


a 
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^2 • (n - C 2 AQ - 02Ag)/u»2 


(7.41) 


0 - (n - A)Aq - YqAq - CqA 2 - Oq 


^ 0^2 ” ***' 


‘oN)' 


A stability calculatloa shows that this solution Is Initially stable, but 
becoaes unstable at a value m ■ > A. 

The solution bifurcating at the point n ■ 4 is given by 


A 2 ■ 0, 


n 


A - 


^ 0^0 



(7.42) 


This is identical with the transcritical solution described in C. of Section 
69 and is represented by the asymffletrlc parabola QPR. However^ the stability 
system shows that the branch PQ is stable (as in Figure 6 )» but that the whole 
branch PR is now unstable. 

In particular we find that there is a value n of wit. 0 < n < A, 

s s 

at which one of the stability exponents for the branch PR changes sign. From 
the corresponding point S on PR there Is secondary bifurcation of a solution 
Joining S to the point T. This Is a mixed-mode solution, and Is actually a 
distinct solution of equations (7.40) < We find Chat this solution Is 
unstable. 

Reviewing the situation as shown In Flgurt; 9, we see that 
before n reaches 0 , the system may stay on the null solution or snap-through 
to Che axisymmetrlc branch PQ; on 0 < n < n^ the system may choose PQ or Che 
predominantly m ^ 2 solution OM, ON, OL; but Chat for n > the only stable 
solution Is Che axisymmetrlc mode represented by the branch PQ. Thus, even 




-64- 


though the m • 2 mode appears first according to linear theory, It soon loses 
out to the axis ymme trie mode as a result of the Interactions. 


We consider next the case where a Is slightly greater than ag, 
a > ag, with 


M - M 


Oil* 


4 - M 


211 


- M 


Oil 


> 0 . 


(7.43) 


The governing equations are <'lth n and n - 4 Interchanged, and the 

stability equations are (7.40) with n and n - 4 Interchanged. The results are 
shown In Figure 10 for large and moderate Prandtl numbers; for very small 
Prandtl numbers the picture Is generally the same apart from reflection In 
the n-axls. 


Figure 10 

The solution bifurcating from n ■ 0 Is the transcrltlcal axlsymmetrlc 
solution given by 


42 - 0, n - TfgAQ - u.qaJ - 0. (7.44) 

represented by the curves PQ, PR In Figure 10. Calculations based on (7.40) 
show that the branch PQ Is stable and PO Is unstable, as In the non- 
Interactlve case depicted in Figure 6, but that OR, Initially stable, loses 
stability at n • where 0 < n <4. This re-emphasizes the preference for 
the upper branch discussed In Section 8. 


The solution bifurcating from n • 1 s a solution of the equation* 
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A 2 ■ (n - A - CjAq - 


0 - hAq - YqAq - CqAj - <>oN)^2 " 0* 


(7.45) 


It la found to exist for n > 4 and to be a slight modification of the pure m ■ 
2 aolutlon deacrlbed In B. of Section 6. This solulon Is represented by the 
curves AL, AM and AN In Figure 10. Calculations show It to be unstable. 

Finally, there Is another solution of (7.45) which Is a mixed mode and 
which bifurcates from the axlsymmetrlc solution (7.44) at the point S. This 
solution Is represented by the curves SU, SV, ST in Figure 10, and Is found to 
be stable. 

We see, therefore, that there are two stable solutions: The pure 

axlsymmetrlc mode PQ, and the mixed mode emanating from S. The latter 
represents a distortion of the lower branch duo to the modal Interaction. 
Nevertheless the upper branch PQ Is preferred as In the non-lnteractlve 


situation. 
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8, DISCUSSION AND CONCLUSIONS 

We have considered Marangonl instability in a circular cylinder under 
the simplifying assumptions that the upper free surface is non -deformable, 
l.e* C and the sidewalls are adiabatic and impenetrable but "slippery". 

The linear stability curves vary with surface Biot number h and Rayleigh 
number R as expected from the analyses of the unbounded layer. (The boundary- 
value problems are identical.) decreases with R and increases with h and 
these features are shown in Figures 3-5. Here, the envelope of each set of 
curves gives for each aspect ratio. Had we used the more realistic 
rigid-side-wall boundary conditions, this envelope would have been modified. 

We would still expect to have interlacing of the modes though the modes might 
Interlace in a different sequence. Only the direct calculation of these 
neutral stability curves could determine this. We shall aiscuss below the 
implications of the use of slippery sidewall conditions. 

Given the qualitative similarities of cases for various values of h 
and R, we Investigate the nonlinear behavior only for the single set h ■ 0, 

R « 0. We have selected five aspect ratios and performed the bifurcation 
analyses for these. 

Cases A, B and C of Section 6 relate to aspect ratios corresponding to 
simple eigenvalues for which m ■ 1,2,0 respectively. We see that for m 0 

that we have supercritical bifurcation only. However, when m * 0, the axi- 
s%Tnnetric convection is subcritical and snap-through convection can be expected. 
It is only in this axi-syninetvlc case that flor direction is distinguished. 

For low Prandtl number Pr, there is downflow in the center while for all 
Pr :: 1 , there is upflow in the center. As is well-known, subcritical Irsta- 
bilitles have associated transport values with hysteresis behavior. (This 
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vlll be discussed further for Case B of Section 7.) It is ussy to calculate 

the degree of subcritlcality possible in the oxl-symnetric mode i.e. the 

tH/H at the nose of the curve. This varies from 0.18% at Pr - 0.1 to 1.4% 
c 

at Pr ■ **. The only comparison available is to the results of Scanlon and Segel 

(1967) who examine an "infintely deep" layer having no sidewalls. In their 

analysis they have Pr > • and find 4 m/M^ " 2.3%. The two analyses are in 

reasonable agreement. This gives further substance to our feeling that our 

results reflect the inherent nonlinear behavior of the system. 

Case A of Section 7 examines a neighborhood of a a^ of Figure 3 where 

M is a double eigenvalue of modes m ■ 1 and m > 2. The nonlinear theory 
c 

gives a coupled pair of nonlinear amplitude equations, (7.15). The analysis 

shows (Figure g) that for a slightly larger than a^ the first mode to appear 

(at M “ M^) is the pure mode m > 2 as predicted by linear theory. As M is 

increased above the system remains in this mode and possibly no further 

transition is predicted. Alternatively, the system may progress through there 

the sequence: pure m > 2, mixed (1,2) mixed time- periodic and perhaps pure 

m « 2 as discussed in Section 7. On the other hand if a Is slightly smaller 

than a. the transition sequence (Figure 7) is completely different. Here, at 
A 

M ■ M , a mixed (1,2) mode occurs and this mode becomes unstable for an M > M . 
c c 

Hence, there must be a transition to the pure mode m ■ 2. We find then that the 

mode m ■ 2 may persist and be stable for M large enough on either side of a « 

independent of the prediction of linear theory. This result depends on the 

stability of the time-periodic mode which has not been examined here but will 

be investigated further in later work. 

Case B of Section 7 examines a neighborhood of a ■ cf Figure 3 where 

M is a double eigenvalue of modes m * 2 and m ■ 0. The nonlinear theory 
c 

gives a coupled pair of nonlinear amplitude equations (7.36). The analysis 
shows (Figure 8) that for a slightly larger than a^ the first mode to appear 
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(ac M ■ M^) is Che pure oxl-syinneCrlc one. There is a snap-through at some 

M<M (Che snap-through is accompanied by a dynamical hysteresis behavior), 
c 

With sufficiently small disturbances the system resides in this m ■ 0 mode 
for all M covered by the theory. Yet, for M large enough there can be a 
transition to the mixed (2,0) mode if the disturbances are large enough. On 
Che other hand if a is slightly smaller than a (Figure 9) very small distur- 

O 

bances evolve as a mixed (2,0) mode until this mode becomes unstable at 
supercritical conditions. At this point there is a snap- through transition 
Co the pure m * 0 mode. If M is then decreased, a dynamical hysteresis loop 
would be revealed since the Jump back to the mixed mode would usually occur 
at much lower values of M; conceivably in fact the jump could be to the state 
of pure conduction. This interesting hysteresis loop could consist of three 
distinct states: mixed (2,0), pure m ■ 0, pure conduction! Alternatively, 

if Che system is very noisy in that large disturbances are present, as M is 
increased from subcricical condition, the system could snap-through directly 
from conduction (at M < M^) to the axi-synnetric state and completely by-pass 
the mixed-mode state for increasing M, yet return to it when M is decreased. 

The above analysis should give a faithful representation of the nonlinear 
processes in fairly small containers in which Marangoni instability cakes 
place. If the replacing of the ideal sidewalls with more realistic rigid 
walls does not change the sequence of the modal inter lacings, then the theory 
could be applied to experiment in a straight-forward way. ^ the. sequence 
of nodal interlacings does change, then the theory should be applied by 
first locating the double eigenvalue by experimental observation. Thus location 
a > a , say, would be different from that of the "slippery" wall linearized 
analysis. However, once it is located, the raising of M for aspect ratios a 
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on olcher tide of should be well represented by the above theory. It is 
thus a relatively "simple" observation of flow pattern that would Initially 
need verification. 

There has been no previous nonlinear analysis of Marangonl convection 
In a bounded container. The present work represents a first exploration of the 
phenomena albeit with an idealized model. The idealization on the upper free 
surface C**0, will be relieved in our future work so that effects of free- 
surface deflection can be assessed. The dropping of idealization of slippery 
sidewalls entails a major computing program that will not be undertaken. 
Clearly, certain small imperfections on either the sidewalls or the free 
surface can lead to an imperfect bifurcation in which the predicted sharp 
instabilities become gradual changes. Our work here provides the framework 
for studying these effects as well. 
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FIGURE 1: 


FIGURE 2: 


FIGURE 3: 


FIGURE 4: 


FIGURE 5: 


FIGURE 6: 


FIGURE 7: 


CAPTIONS FOR FIGURES 

A sketch of a possible static llquld-gas configuration in a wide 
container (a) on Earth at Ig, (b) in space at 10"^g, (c) in space 
at 10 ^g, (d) in space at 10”^g, 

A sketch of a possible static liquid-gas configuration in a 
narrow container (a) when the interface is curved and (b) when 
the interface is flat. 

Calculated stability curves, M versus a, for h ■ 0, R ■ 0 where 
m and i are the azimuthal and radial wave numbers, respectively. 

is the envelope (minima) of the curves. 

Calculated neutral stability curves. versus a for h ■ 0 and 
various values of R. The m values take the same sequence as in 
Figure 3. 

Calculated netural stability curves, versus a for R * 0 and 
various values of h. The m values take the same sequence as in 
Figure 3. 

The bifurcation diagrams for a - 1.80 corresponding to a pure 
axl -symmetric mode m >■ 0. Solid lines represent stable branches 
while dotted lines represent unstable branches, (a) The case 
of Pr ■ 0.1" the branch PQ represents downflow in the center. 

(b) The case of Pr 2 1; the branch PQ represents up flow in the 
center. 

The bifurcation diagrams for a slightly less than a^ s: 1.20 where 
^ ■ M 2 ■* M^. Solid lines represent stable branches while dotted 
lines represent unstable branches. 


-72 


riCUIE 8: 


FIGURE 9: 


FIGURE 10: 


The blfurcaclon diagrei.e for e sllghcly grMCer than 1*20 

where A - • M^. Solid lines represent stable branches while 

dotted lines represent unstable branches. The curly lines represent 
tlae-periodlc bifurcstions. 

The bifurcation diagrams for a slightly less than a^ sr 1.70 where 
A ■ - Mq. Solid lines represent stable branches while dotted 

lines represent unstable branches. 

The bifurcation diagrams for a slightly greater than a^ 1.70 
where A ■ Mq - M 2 . Solid lines represent stable branches while 
dotted lines represent unstable branches. 







/ / 
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FIGURE 7 
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CHAPTER II . CONVECTIOM IN RECTANGULAR CYUNDERS 

1. INTRODUCTION 

In Chapter I, which we shall call I, we discussed Marangoni instabilities 
in a circular cylinder and distinguished between simple eigenvalues and 
double eigenvalues, where secondary bifurcations are possible. 

In the present chapter, we examine Marangoni instability in rectangular 
containers. Again, we assume that the upper free surface is non-deformable, 
C^O, and Che sidewalls are adiabatic and impermeable but "slippery", which 
in Che rectangular geometry, corresponds Co zero shear stress applied on Che 
boundary. We use Che asymptotic theory of Rosenblat (1979) Co examine the 
steady convective states near and how transitions from one state Co another 
occur. We limit ourselves Cc interactions of modes in the form of two- 
dimensional roll -cel Is, which are predicted for rectangular containers having 
Che shorter side comparable to Che depth and Che longer side larger than Che 
depth. Since much of the full development is similar to that in 1, we give 
only those details which distinguish convection in a rectangular container 
from convection in a circular cylinder. 
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2. FORMULATION 

Consider a viscous liquid, which partially fills a container 
of rectangular cross-section. The mean depth of the liquid is d and a 
horizontal cross-section has lengths a^d and a 2 d in the x and y directions 
respectively. Hence, a^ and a^ are the aspect ratios. The axis of the 
cylinder is anti-parallel to the direction of gravity, and the upper surface 
of the liquid is open to an ambient gas. 

The development of the non-dimensional nonlinear disturbance equations 
and boundary conditions parallels that in I. Again, in the limit of 
small capillary number and when the lateral boundaries are adiabatic 
and impenetrable but "slippery", we obtain the following nonlinear problem. 
From equations (2.19-2.21) of I, 

Pr Hl-f ^ + (v*^)v|’ “ -7p + V^v + M (2.1) 

7«v “ 0 (2.2) 

M ■'''II - w + (v7)e} - v^e (2.3) 

where M, R and Pr are the Marangoni, Rayleigh and Prandtl numbers defined 
in equations (2.22) of I. 

The bottom of the box is a rigid perfect conductor so that 

0*0 on z*0 , , OsySa, , (2.4a) 

v*0onz*0 , , 0sysa2 . (2.4b) 

Since the capillary number is zero^the upper surface is flat, so that 
from equations (3.5) of I, the heat transfer condition is 
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(2.5a) 


|^ + h0-O on z-1 , Osxsa^ , 

and the kinematic condition is 

w*0 on z*l t Oixsa^^ , 0^ySa2 . (2.5b) 

The conditions of thermocapillarity become 

u+w+e - v+ w+0 »0 , z«l , Oacx^a. , O^y^a. 
zxxzyy i ^ 

(2.5c) 

Finally, the "slippery" sidewalls reduce in Cartesian coordinates to 
adiabatic, impermeable stress-free planes. These conditions take the form 


u»w "V *0 *0 on 

XXX 

9 < 

1 

o 

0 :: y :: a^ 

, 0 i z i i 

(2.6a) 

v*w *u *0 '^0 on 

y ■ 0 ,a 2 *, 

0 • X a^^ 

, 0 < z i 1 

(2.6b) 

y y y 



The critical Marangoni number at which the conduction solution loses 


stability is determined from linearization of the syztex (2.1)- (2.3) 
together with the (linear) boundary conditions (2.4)- (2.6). As in I, we 
assume that occurs at a state of neutral stability so that the governing 
linearized equations bcctxne 

V^w + - 0 (3 

and 

7^9 + Mw » 0 . (3 

System (3.1) plus boundary c'>nditions (2,4)-(2.6) may be solved by 
seeking separable solutions in the form 

w(x,y,z) » cos C mj^TTx/a^] cos [m 2 ny/a,JY(z) (3 

0(x,y,z) « cos L m-j^nx/aj^] cos Cm 2 ny/a 2 jX(z) (3 

with similar definitions for u and v. Here X and Y are the same functions 
as those in I since they are eigenfunctions of a set associated with the 
stability of the infinite conductive layer considered by Pearson (1938). 
Here m^ and m^ run over all non-negative integers. 

When forms (3.2) are substituted into eqns. (3.1), an effective wave 
number X appears where 

■ [(fflj^/a^)" + (m2/a^)^]rr^ . (3 

The effects of buoyancy through the Rayleigh number R and the effects 
of the free surface being a poor insulator through the surface Biot number 
can be explored as in I. The effects are the s.une in that increasing R 
decreases M and increasing h increases M . Tlicse results will not be 





We note that for an infinite layer, X is the overall wave number t^ich takes 
on all values on [0,“), M(X) would then have the minimum * 79.6 for 
X 2. This result is due to Pearson (1958). For the present enclosed 
layer, M(X) must be minimized over only those admissible X given by eqn. (3.3). 

The relationship between the box aspect ratios and the mode of 
convection, indicated by the integers (m^tm^), is given implicitly 
by equation (3.4). We have evaluatec this relationship for a range of box 
sizes for all possible modes of convection. The results are given in 
Figures 1-6, in which M is given as a function of for fixed values of 
For clarity, modes having large critical Marangoni numbers are not 
shown. Consider firi<>t the case of a^ * 0.5 shown in Figure 1. As the box 
size, a^ , increases, the preferred mode, i.e. the mode having the lowest critical 
Marangoni number, changes in a specific way. This sequence is among modes 
for which m 2 * 0. Thus, we have two-dimensional roll cells whose axes are 
aligned with the shorter dimension of the box. We shall call these x-rolls. 

It is seen that for box sizes a^^ ■ — ^ with mj^ » 1,2,3»»* that X ^ 2 and 
the critical Marangoni number is minimum at the value 79.6 appropriate 

to infinite layers. Away from these values, the sidewalls exert a 
stabilizing influence, even though they are "slippery". While the fact that 
several box sizes can have the same M * is presumably an artifact of 
the use of the slip-wall boundary conditions, the existence and progression 
of preferred modes due to the finite size of the container is not. 
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In the case of buoyancy-driven convection, the dependence of M on box size 
is monotonlc but has kinks at the points of moda switching, Davis (1967), 
and the effect of sidewalls is Co align the roll axes with Che shorter side 
of Che box. This is the same progression and alignment as predicted here 
for ^2 * shall see below that the present treatment leads to 

some predictions of preferred mode orientation which are presumably arti- 
facts of the slip-wall boundary conditions. 

To sumnarize, the curves in Figure 1 predict preferred modes consisting 
of x-rolls, and Che progression is to add more x- rolls as the box size 
increases. Of particular interest are the aspect ratios at which two modes 
have Che same critical M, this is a double eigenvalue of the linear theory. 

Figure 2 shows Che results for a^ ■■ 1.0. Since the modes with * 0 are 
unaffected by the length a^ the lower curves are identical to chose of 
Figure 1. We anticipate however, that ae d 2 approaches there are 

two-dimensional rolls with axes aligned with the longer side of the box, 
(y-rolls), which might have lower critical Marangoni numbers Chan Che x-rolls. 
This is not yet the case for the (0,1) mode for the conditions of Figure 2, 
but becomes so for the conditions of Figure 3. Finally, we note the occurrence 
of more complex three-dimensional modes of convection, e.g. the (1,1) and (2jl) 
modes, having Marangoni numbers close to, but above chose for x-rolls. 

Figure 3 gives results for a^ * 1.5, and shows several complex features. 
First we note that the (0,1) y-roll has Ma ^ Ma^ for this value of 
independent of a^, and hence is often the preferred mode. However, since 
a 2 t r/2, there are small ranges of box sizes located near a^ ■ nmj^/2, for 
which x-rolls are preferred. We also note that three-dimensional mode:, 
e.g. the (1,1) and (2,1), become closer to being preferred. At a 2 ■ 2.0, 
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I 


Che results shown in Figure 4 indicate that the y-rolls (0,1), (0,2), r.re 
no longer preferred, and the three-dimensional (1,1) and (2,1) modes arc 
preferred over x-rolls for some range of values of a^ away from a^ > m^Tr/2. 

For 32 * 3.0 (l.e. close to n) Figure 5 shows that a situation analogous 
to that In Figure 3 occurs; the y-roll (0,2) has Ma as Ma^ and is preferred 
for all box sizes a^ away frem m^n/2. Finally, as shown in Figure 6, as 32 
increases, the number of modes competing and having Ma > Ma^ increases, 
and the envelope of these neutral curves becomes nearly the horiconCal 
line Ma ■ This reflects the diminished effect of the sidewalls in 

determining the preferred mode. 

The results may be suimarized by a map in the 
preferred modes. We note Chat the pattern of preferred modes must be 
anti-synmecric about a^^ » a 2 » corresponding to a rotation of Che coordinate 
system. Thus, M(a^,a.,) M(a 2 »aj^) and the preferred modes, (mj^(aj^,a 2 ) , 

m 2 (aj^,a 2 )) “ (m 2 (a 2 »aj^), mj^(a 2 ,aj^)) . It is clear from the previous discussion 
that this map will be complex, and that as a^ and a^ become large, man} 
modes will have values of the critical Marangoni number close to Chat for t 

the preferred modes. This map is shown in Figure 7. With one exception, 
it is difficult Co speculate on the degree to which this complexity depends 
upon the use of slip-wall boundary conditions. Complexity of this degree 
docs not occur for buoyancy-driven convection in a box, Davis (1967), but 
does occur for buoyancy-driven convection in a bounded porous media; Beck 
(1972). The persistence of y-rolls and x-rolls at a^ >« mj^Tr/2, a^ « 
m^ ,m. ■ 0,1,2, ••• respectively, v.i.11 not occur if more realistic no-slip 
boundary conditions are applied. Careful study of results similar to 
those in Figures 1-6 indicates chat much of the co.ttplex mode-switching 
is due to neutral curve for y-rolls (x-rolls), oeing a horizontal line. 



P; 


[ 
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intersectlng many cimes the neutral curve for other modes. No-sllp sidewall 
conditions do not admit pure x-rolls or y-rolls, with the result that Che 
neutral curves for modes whi«.h are close to y-rolls, (x-rolls), may not 
exhibit as many intersections. However, this does not imply that Che 
bifurcation theory developed below will be simpler necessarily, as these 
modes may continue to be near-^nulciple eigenvalues of Che linear theory. 


I 


1 


c 
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4. EIGENFUNCTION EXPANSIONS 


In the nonlinear theory we focus on certain special interactions 


appropriate to one horizontal box-diraunsion being comparable to the depth 


and the other much larger. In particular wc shall take a^ * 1.0 so that 


only x-rolls are predicted by linear theory. It is the interaction of rolls 


that we shall address. Although we must develop the theory for Rayleigh 


number R 0 for completeness properties of the differential system, we 


shall, with no loss of generality, set R » 0 at the end. Hence, pu.e 


Marangoni instability will be analyzed. 


Let us restate the linear stability problem for the case at hand: 


2 -1 
7 V - 7p + M R0^ = 0 


Vv = 0 


7“0 + ^!w = 0 


9 = u = w = 0 


on z = 0 


0 »w*u+9 =0 

Z Z X 


on z = 1 


6 * u * w =0 

X X 


on X = 0,a, 


The problem is now a two-dimensional problem since we are interacting only 


x-.olls; hence, v, ^ = 0. 


For fixed M the eigenvalues are denoted by with m,j » 1,2,*** and 


m is the horizontal wave number while j is the vertical wave number. 


J 


Define 
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The elgenfuncclons are 


7‘V *0 


* ★ 

where X , Y satisfy system (5.14)-(5.16) of I with h 


X ■ m/a. . 

m I 


(4.2) 


u , ■ - X sin X X DY ,(r) 
mj m m mj 


w , ■ cos X X Y .(z) 
mJ ro mJ 


®.J • V 


(4.3a) 

(4.3b) 

(4.3c) 


where the and Y^^ are the elgensolutlons of the system 
(4.10)-(4.12) of I when M - M, R - R ,, X - X and h - 0. 

mj ID 

The adjoint problem Is 

^ it 

7“v - 7p + MO^ ■ 0 


(4.4a) 


2 * »-l * 

7 6 + M Rw - 0 


(4.4c) 


★ 


★ 





9 


u 

■ w 

- 0 

on 

z 



Hr 

He 

★ 



0 

■f 

V 

■ W 

- u - 0 

on 

z 

- 


Z 


z 



★ 


★ 

* 




u 

■ 

V 

« 0 

0 

1 

on 

X 



X 

X 





(4.4d) 

(4.4e) 

(4.4f) 


The adjoint eigenfunctions arc 

u* - - X sin X X DY*,(z) 
mJ m m mj 


w*. • X^cos X X Y (z) 
mj m m mj 

* * 

0 , • cos X X X ,(z) 
mj in mj 


(4.5a) 

(4.5b) 

(4.5c) 
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We now decompose all dependent variables into horizontal mean (i.e. 
x-mean) plus departures from the mean as follows 

v«v + v' t 9*0+6* , p"p+p' (4.6a) 

where for each quantity g, 

_ 1 p‘l 

g • “ g dx . (4,6b) 

®1 **0 


For the case R ■ 0, the equations (2.1)* (2.3) are 
2 -1 

7 V - Vp “ MP^ [v^ + (v*7)v} 

7»v - 0 

7^0 + Mw » M{0j. + (v*7)9} . 


(4.7a) 

(4.7b) 

(4.7c) 


If forms (4.6) are introduced into system (4.7), then we get 

7^ “ Vp = (4,8a) 

7»v » 0 (4.8b) 

79+Mw»m|9^ + (v*7)9 + (v'*7)9'| (4.8c) 

and 

v^v' - 7p' ■ (v'‘V)v + (v*7)v' + [(v'*7)v']^| (4.9r.) 


7*v' *0 (4.9b) 

7“9 ' + Mw' »m| 9^+ (v'*7)0 + (v*7)0' + [ (v' *7)9 '] ^[ (4.9c) 

where [ denotes the fluctuating part of [ I, The same boundary conditions 
hold for both systems (4.8) and (4.9). 

As is well-known, there is no mean velocity field induced by the 


convection and thus 
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V - 0 . (4.10a) 

Equation (4.8c) Chen simplifies considerably, and for steady or quasi- 
Static convection I 

0^ - M(w'e') . (4.10b) 

If these relations are used Co simplify (4.7) we obtain. 


2 

V V* - Vp ' 


MP”S'v^ + C(v'*7)v']jj 


7*v' * 0 


7^9 ' + Mw' 


m( 0' + Mw'(w'0') + C (v'*7)0 '].} 

V t ^ • 


(4.11a) 

(4.11b) 

(4.11c) 


We now take the scaler product of (4.lla,c) with th*. adjoint vectors 

if if 

(^j,8^^j) at M » M^, R * and integrate over the fluid volume. This gives 


(M 


-M Xe^.w') - m"^r .<w",0*> = M(e" 0' + p'V ‘v') 
c ' mj c mj mj mj t r -mj -t 


-1 * 


+ M<0*.i!!(v'* 7)0'],- + Mw'iw'O’)!- 

mj V — I ^ 


+ P ^v , •[ (v‘ *7)v'] -) (4.12) 

r -mj ^ — ■'f 

for each m and j. Equation (4.12) is the basis for the derivation of the 
amplitude equations. 



J 
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5. SIMPLE AND DOUBLE INTERACTIONS 
A. Simple eigenvalue for (1,0) « Let us consider an aspect ratio 
a^^ « 1.5 which corresponds in Figure 2 to a simple eigenvalue for 

convection with (mj^,m 2 ) ■ (1|0). 

We find that 


M - 79.4 
c 


(5.1) 


and by hypothesis, « 0. The quadratic interaction of mode 11 get crates 
the 21 mode (R 2 X ^ set.^is 


. {11,21}. 


(5.2a) 


We write 


(0’,v') * ^ 2 ^® 21 ’-21^ 

and substitute into equation (4.12). We obtain 


(5.2b) 


''1*1 ■ - h 

“2^ - - - "2 


when M = M , where 
c* 


V = d"Hl <9* 9 ,+Pr”^v*. *v > 
m m c ml ml “ml -ml 


f = d ^<w ,9 , ) 
m m ml ml 


(5.3a) 

(5.3b) 

(5.3c) 

(5.3d) 


d * (8 , w , 
m ml ml 


(5.3e) 


We shall not give all the details of the evaluations since they are parallel 
to those of I. After a good deal of algebraic manipulation, we find that 
if a, and are constants that 


Zi - ttj^Aj^A, + Bj^Aj^ 


(5.3f) 
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^2 - 


and the governing amplitude equation takes the form 




The computations of the coefficients have been performed for various values 
of Prandtl numbers and some results arc shown in Table 1, 


Pr 

-4 

X 10 

X lO” 

0.1 

0.37 

5.2 

1.0 

0.13 

0.64 

10.0 

0.10 

0.43 

a 

0.10 

0.41 


Table 1 


Since > 0, the convective state results from supercritical bifurcation 
and is stable. 


B. Simple eigenvalue for (2,0) . Let us consider an aspect ratio 
* 3.1 which corresponds in Figure 2 to a simple eigenvalue for 

convection with (ra. ,m^) ^ (2, OK We find that 


M -79.2. 
c 


There is again a quadratic interaction anu the sctxf is 
sf - (21,40 

We ODiit all details and scati. the final amplitude equation, 


• (M - M J.V 


<) “ ^*>^7 


J 




where the coefficients have the numerical values given in Table 2, 


Pr 

V 2 X 10 ”^ 

u)2 X 10 

0.1 

0.44 

1.4 

1.0 

0.13 

0.26 

10.0 

0.10 

0.25 

OD 

0.098 

0.25 


Table 2 

Again, > 0, the convective state results from supercritical bifurcation 
and is stable. 

C. Double eigenvalues for (1.01 and ('2. OK Let us consider an aspect 
ratio aj^ ■ 2.21 which corresponds in Figure 2 to a double eigenvalue for 
convection with * (1,0) and (ni^,m 2 ) “ (2,0). We find that 


M » 90.2 
c 


(5.8) 


The quadratic interaction of modes 11 and 21 generate modes 31 and 41. The 
sety<$ is 


- (11,21,31,411. 

We write 

W.v’) 

1*1 

and substitute into equation (4.12). We obtain 
• (M-M^)A^ - 

v 2^2 “ ■ ^2 

^^*^ 3 ^ 3 ^ “ • "3 


t5.9a) 


(5.9b) 


(5.10a) 

(5.10b) 


(5.10c) 
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”c\V4 “ ■ h <5-lW) 

wheu M ■ M . Here Che £., f , , R., R. and functionals Z, - Z. are defined 

c 3 4 3 4 14 

in analogous way as in J. Again, we ooic details and state the final 
amplitude equations; 


A ■ ■ ^*1*1 

(5.11a) 

2^ ■ W-M..)*2 - <4*? - ®2*1*2 ■ “ 2 '^ • 

(5.11b) 


Numerical values of the coefficients are given in Table 3. 


Pr 

X 10 "^ V 

2 X 10 "^ 

-3 

X 10 

02 X 10 "^ 

3^ XlO" 


0 ^ X 10'^ 

X 10‘^ 

0.1 

0.58 

3.9 

0.88 

1.5 

0.18 

2.3 

1.3 

5.0 

1.0 

0.23 

1.2 

0.34 

0.43 

0.18 

0.33 

0.17 

0.95 

10.0 

0.20 

0.90 

0.29 

0.32 

0.18 

0.24 

0.12 

0.77 

OD 

0.20 

0.87 

0.29 

0.31 

0.18 

0.23 

0.12 

0.76 

We analyze the 

equations 

(5.11) 

Table 
in detail 

3 

below. 





D. Double 

eigenvalue for 

(2,0) and 

(3,0). 

Let ua consider an 

aspect 


ratio a^ > 3.81 which corresponds in Figure 2 to a double eigenvalue for 
convection with (m^^.m^) - (2,0) and (mj^.m^) “(3^0)^ We find that 


M » 82.9 
c 


(5.12) 


The quadratic interaction of modes 21 and 31 generates modes 11, 41, 51, 61 
so the sec^ is 

a1 - [11,21,31,41,51,61]. 


(5.13a) 


We vrlte 


(e',v') - j A <9 ) 

and substitute into equation (A. 2). We obtain 
V2^ - (M-ri^> - 


V3 * • h 




when M • Here the R^, Z^ are defined in analogous way as in I. 

Rather than give the details, we state the final .unplitude equations. 

V2A2 - (M-M^)A2 - u»2A^ - T2A2A5 
V3A3 - (M-M^)A3 - T3A2A3 - 
and the coefficients aie given in Table A. 


Pr 

X 10 ’^ 

V3 X lO”^ 

* A 

ui2 X 10 

T2 X 10‘^ 

0)3 X 10‘^ 

T3 X 10 

0.1 

0.A6 

3.6 

l.A 

- 7.8 

5.1 

6.7 

1.0 

0.18 

1.2 

0.17 

- l.l 

0.85 

1.2 

0.0 

0.15 

0.92 

0.12 

- 0.76 

0.66 

0.89 

OB 

O.IA 

0.89 

O.ll 

- 0.72 

0.6A 

0.86 


(5.13b) 

(S.lAa) 

(S.lAb) 

(5.1Ac) 

(5.15a) 

(5.15b) 


Fable A 
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6. ANALYSIS AND DISCUSSION 

In cases A and B of Section 5, the self-interaction of roll cells (1,0) 
and (2,0) is considered. In both cases, the interaction is governed by 
single amplitude equations containing cubic but no quadratic noiilinearities. 

These are equations (5,4) ami (5,7) respectively. The values and V 2 , 

depending on Prandtl number Fr, are values from the linear stability problem 
and for given X of equation (3.3) are identical here to those of I. Careful 
comparison shows this. The values and {^2 of equations (5 4) and (5.7) 
are always positive so that these simple self-interactions always correspond 
to stable supercritical bifurcation. It is easy to show that for any values 
(m^^.m^) ^ (0,0) indicated in Figure 7 that self-interactions always have 
amplitude equations of the sar.ie form, i.c. 

vA » (M-M 1 a - uuA (6.1) 

c 

where v > 0. Presumably, u) > 0 for any of these so that stable, supercritical 
bifurcation is always predicted for self-interactions. This is likewise 
true in the case of the circular container of I for m ^ 0. It is only for 
the (m ■ 0) axisynroetrlc mode that equation (6.1) is augmented by a quadratic 
term. Thus, the axisynmetrlc mode blfurcctes subcritlcally and so has 
snap-through and hysteresis properties as discussed in I. 

In case C of Section 5, the interaction of modes (1,0) and (2,0) is 
exaalned near the double eigenvalue at a^ - 3.1 of Figure 2. The tjoveming 
amplitude equations (5.11) .ire a pair of coupied equations identical in 
form to equations (7.10) md (7.11) of 1 which govern the interaction if 
modes m ■ 1 and m ■ 2 rear their double eigenvalue. Again, the are 
linear theory values th-t depend on i’r and \ but not on the cylinder gecro try. 

J 
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A1 though Che coefficients are not identical In Ch-t two cases, all of the 
qualitative predictions are. Figure 8 shows Che results of our analysis 
of (7.10)-(7.11) for a^ < 3.1. The mixed mode containing both modes (1,0) 
and (2,0) bifurcates supcrcritically at and as M is further increased, 
fellows either or ^2 *^2^2* * value of 

T| ■ greater than A a M 2 “Mj, there is secondary bifurcation to a pure 

mode m * 2. This branch is labelled T^S. 

Figure 9 shows the situation for > 3.1. Here, as M crosses M^, Che 
pure mode m * 2 bifurcates supercritically and follows either curve 
OjT^ or O^S. However, for T) a M - less than A a pure mode 

persists but only on the branch O^S. Again, there is the possibility of 
branch O^T^ bifurcating first to the mixed mode and Chen to time periodic 
convection. The amplitude equations (5.11) arc in form identical to Chose 
governing hexagonal cells as predicted by Scanlon and Segel (1967) for 
horizontally unbounded layers. However, since the contexts are quire dif- 
ferent, the coefficients are quite different. Scanlon and Segel find 
subcritical hexagons. We find only supercritical convection of mixed mode 
or pure mode m ■ 2. 

In case D of Section 5, the interaction of modes (2,0) and (3,0) is 
examined near the double eigenv/alue at ■ 3.81 of Figure 2. The governing 
amplitude equations (5.14) are a pair of coupled equations. Again, are 
linear theory values that depend on Pr and X but not oti the cylinder geometry. 
Figure 10 shows the situation for a^ < 3.81. 'Ihe pure mode (2,0) bifurcates 
supercritically at M^, T] a ■ 0, and steady convection follows either 

branch OS as M increases. At a value of ■" > A a ^2 ” *^l secondary 

bifurcation to the mixed mode containing both modes (2,0) and (3,0) and as 
H increases further, A 2 follows either branch SU and A^ follows a branch ST. 
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Figure ll shove the situation for > 3.81. Here, at the pure mode 

(3,0) bifurcates supercritically and follows either branch OS until 
T) ■ M - < A ■ “Mj. Here, there is secondary bifurcation to a 

nixed node in which A2 follows either SU and follows an ST. The sequence 

of events here, near a ■ 3,81, has no counterpart in I since there was no 
double eigenvalue there for nodes m ■ 2 and m > 3. However, the amplitude 
equations (5.14) have the form typical of Rayleigh-B/nard convection in 
containers as discussed by Rosenblat (1981). Hie Figures 10 and 11 are 
typical of Rosenblat 's results idilch apply to the buoyancy-driven convection. 

In sunnary, we again find that interactions near double eigenvalues 
give qualitative features that strongly distinguish behavior for aspect ratios 
on one side from behavior on the other side. Parallels as well as diffcrenc-i 
in behavior exist between the circular and rectangular cases. 


t 
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CAPTIONS FOR nCURES 

I: Scabllicy curves M versus for L > 0 ac S 2 ■ 0.5. The 

pairs ( 0 ^, 02 ) denote integral number of cycles in (a^,a 2 ). 

2: Stability curves M versus a^ for L - 0 at a 2 ■ 1.0. The 

pairs (m^.ffl 2 ) denote integral number of cycles in (a^,S 2 ). 

3: Stability curves M versus a^ for L * 0 at a 2 ■ 1.5. The 

pairs (m^,m^) denote integral number of cycles in ( 0 ^, 02 ). 

4: Stability curves M versus a^ for L ■ 0 at a 2 * 2.0. The 

pairs (a^(ffl 2 ) denote integral number of cycles in (a^,S 2 ). 

5: Stability curves M versus a^ for L ■ 0 at a 2 • 3.0. The 

pairs (m^,m 2 ) denote integral number of cycles in ( 0 ^, 02 ). 

6 : Stability curves M versus a^ for L > 0 at 82 ■ 3.5. The 

pairs (m^,m 2 ) denote integral number of cycles in ( 8 ^, 02 ). 

7: Stability map for preferred mode as a function of a^ and 82 . L = 0. 

The figure is symmetric about - 82 . The pairs (m^. 02 ) 
denote Integral number of cycles in ( 0 ^, 82 ). 

8 : The bifurcation diagram for 82 ■ 1.5, and a^ slightly less than 

3.1 where A « M 2 - Solid lines represent stable branches 

vdiile dotted lines represent unstable branches. 

9: The bifurcation diagrams for a^ • L.5, and a. slightly greater 

than 3*1 where ^ Solid lines represent stable 

branches while dotted lines represent unstable branches. 

The curly lines represent time-poriodic bifurcations. 

10: The rifurcati^Q diagrams for * 1.5, and a^ slightly less chan 

J.8I where ^ - M 2 . Solid lines represent stable branches 

while dotted lines represent unstable branches. 

11: The bifurcation diagrams for a., ■ 1.5, and slightly greater 

than 3.81 where Solid lines represent stable 

branches while dotted lines represent unstable branches. 
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APPEND} /: A 






The classical Raylelgh-Benard problem for Che stability of a 
horizontal layer of Bousslnesq fluid heated from below can be reduced to Che 
solution of the pair of ordinary differential equations^ 

.2 2 

^ - a^)6 + w ■ 0, ” a^j^w - Ra^9 *0, (1) 

dz dr 

where 8, w are respectively Che temperature perturbation and the vertical 
velocity-component perturbation; a Is Che wave number of disturbances In the 
horizontal plane; z Is the vertical coordinate; and R Is the Rayleigh number, 
proportional to Che vertical temperature gradient In the equilibrium state* 
Equations (1) hold on the Interval 0 < z < 1, and their solutions are required 
to satisfy prescribed boundary conditions of Che form 

P^6 + Q^w *0 on z ■ 0, P^0 + Q^w *0 on z ■ I, k ■ 1, 2, 3. (2) 

Here P^, Q^, are linear differential operators whose particular forms 

depend on the characteristics of the boundary surfaces* 

For a layer of unbounded horizontal extent Che buundary-value problem 
(l)-(2) has u nontrivial solution for infinitely many values of R* More 
precisely, there Is a countably Infinite number of simple eigenvalues for each 
wave-number a given by a characteristic equation of Che form 


with 


R • R (a), n • 1, 2, 3, . ** 
n 


Rj(a) < R^fa) < Rj(a) < *** 


(3) 

(4) 





1 






i 
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for each a« The critical Rayleigh nunber for the onset of Instability la 
given by 

R ■ R,(a ) ■ min R (a). (5) 

c 1 c 1 

a 

Corresponding to each R^ there Is an elgensolutlon 

0 (z; a), w (z; a), n - I, 2, 3, (6) 


I 


for various boundary conditions of the form (2) It can be shown that the 
eigenvectors (6) constitute a complete set, a fact which may be useful for 
computations of the corresponding nonlinear stability problem* 

As an example, the simplest case Is that when both boundaries are 
Isothermal and stress-free. Conditions (2) then become 



d^w 

dz^ 


- 0 


on 


0 and z ■ I, 


(7) 


whereupon the characteristic equation (3) has the well-known form 


^ [(mr)^-Ha^]^ 

^n 2 

a 

with eigenvectors 


( 8 ) 


6 ■ sin miz, w ■ [(nn) + a Jsln nnz. (9) 

n n 

o 

Pearson- studied the Marangonl problem, In which the Instability was 
driven by surface-t >nslon gradients rather than by the buoyancy force. 
Subsequently, Nleld^ Investigated the problem allowing the presence of both 




i 
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buoyancy and aurface tension. In these cases the governing equations again 
reduce to the system (1) (with R ■ 0 in Pearson's work). The boundary 
conditions have the general form (2) but of necessity contain a new 
paraoieter» the Marangoni number M. In the specific problem considered by 
both Pearson and Nleld the boundary conditions were 


and 


dz 


3w - 

— ■ 0 on z 
3z 


d w , ^ ^ 

— r + Ma 0 ■ 0 on 

dz^ 


z - I, 


( 10 ) 


(ID 


where h is a surface Nusselt number. These correspond to a rigid, isothermal 
lower boundary and a stress-’f ree, conducting upper boundary. 

Both Pearson and Nield were concerned to determine a critical 
Marangoni number for the onset of instability for fixed values of the other 
parameters of the problem. Detailed calculations by Nield^ led to a single-* 
valued characteristic equation of the form 


M - M(a, R, h); 


( 12 ) 


in other words, the boundary-value problem (1), (lO)-(ll) has a unique 
eigenvalue M for each value of a, R and L (and, correspondingly, a unique 
eigenvector). The critical Marangoni number was then found to be 


M - M (R, h) - M(a , R, h) • min M(a, R, h) 
c c c 


a 


(13) 
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wlth ■ a^(Rf h)* example, Nleld^ found M^(0, 0) • 79.607 with 

•^.(0, 0) - 1.993. The graph of the function (12) for the range 
-500 < R < 1000 with a ■ 2 and h ■ 0 la shown In Figure 1. 


Figure I 


There appears to be t discrepancy between the Ray lelgh-Benard problem 
and the Marangonl problem: the former has countably many eigenvalues while 

the latter has only one eigenvalue M. ‘^is Is somewhat surprising, since the 
boundary-value problem (1), (lO)-(ll) incorporates both buoyancy and surface- 
tension effects. A possibly disturbing consequence Is that the Marangonl 
problem does not seem to have a complete set of elgenfunct lons^ which could be 
usud for computations of nonlinear stability. 

This discrepancy Is reconciled by noting that the Inverse of the 
function (12) has countably many values of R for each fixed value of M, a and 
h. In fact, the Inverse of (12) is of the form 

R - R (a, M, h ), n - 1. 2. 3. ... (14) 

n 

Thus, for each M there are Infinitely many eigenvalues R^^, but only cue 
eigenvalue M for each R. 

We have computed the functions (14) by Inverting the explicit formula 


(12) given by Nleld. In our notation 
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M - A/B 


(15) 


where 


- - d^-R «• (nm)^d .2 - 

2 r m 1 r ' m ^ 2Rd*^[ 'I l™Li_Js2SSL j ^ (16) 


J-J. I ^ 

m* I m m" 1 m 


m*l m 


2 2 

^ , v2 • (nrr ) d ^ / \2 « (ran) d cosmn 

2a^ ^ ^ V (ran) cosrair v ra 


1 A A 

> I ra m* I m 


-2.^ ): 


ra*l m 


I 


(17) 


tn“l 


with 


2 2 3 2 

d • (am ) +a, A “J -Ra. 

m mm 


(18) 


CaloulaClona were performed for the case h ■* 0 and a ■ 2, and are shown In 
Figure 2 as a graph of R against M, different scales being used because of the 
large range of parmeter values. The dotted segment of the lowest curve 
corresponds to the curve depicted in Figure 1. 


Figure 2 


We find from the calculations that the curves Intersect the R-axls et 
values Rj - 67b, R 2 ■ 2.11 >« 10^ and R 3 - 2.05 >« 10^. These compare well with 
the values obtained from formula ( 8 ) at a * 2 , namely, R^ * 668 , 

R 2 • 2.06 X 10^ and Rj - 2.00 * 10^. Also, the lines R - 6.7 “ 10^ and 
R • 9.9 ^ 10^ are asymptotes for the curves. 
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This analysis shows that although a critical value of M can be 
determined, the Marangoni number is not strictly speaking an eigenvalue of the 
boundary-value problem, whereas the Rayleigh number is. 
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CAPTIOWS FOR FIGURES 

The first critical curve R versus M for h 0. 

The first three critical curve R versus M for h ■ 0. 
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APPENDIX B 

The eigenfunctions used in the expansion procedure are 

the solutions of the boundary -value problem (4. 10) -(4. 12) for fixed M ■ 
for R * R . . and with h * 0. Using the Fourier series method of Nield (1964) 

ml j 

we find Chat 


X . .(z) ■ 2 + y a sin nnz 
mil 1 n 

n«l 




nnz 


n-1 


( 1 ) 


where Che coefficients a^, are defined in the following way. Let 
"n ■ 

A - d^ - R . 
n n mij mi 


( 2 ) 

H) 


and 


J - d^ - R , , - (nr)^M 
n n miJ c 


W 


Then 


2 2 
2X .cos ITT J - (nn) Q 
ml n 

irrA 


(5) 


and 


3 - 

n 


nrr a a - 2X .cos nn 
n n mi 

nrr M 


( 6 ) 


where 


1 2X 


mi 




k-1 


^ OP 

XI'* (krr)^cos krr/A, } 


k-l 


( 7 ) 
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The formulas (Al) can also be used co calculate the first derivative of X . , 

m j 

and the first and second derivatives of wh'.ch are required in the nonlinear 

interactions of Sections 5-7. 

* * 

The adjoint eigenfunctions X^^^ and are the solutions of the adjoint 

boundary-value problem (5. 14)- (5. 16) with M * R ■ h ■ 0, as before. 

We find that 


Q0 U0 

X*.,(z) ■ z + y Y sinnnz , Y (z) - V 6 sini 
mij n mij n 


where 


2\^.|(d^ - R . ,)cos nn + (nn)^R . .hI 
mil n mil -■ n mil , 


2nnM (d H - cos nrr) 
c n 


OD 

^ ^ (kTr)^cos 




